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ABSTRACT 

We use infrared spectroscopy and photometry to empirically define the intrinsic, thermal in- 
frared spectral energy distribution (i.e., 6-100 [a.m SED) of typical active galactic nuclei (i.e., 
2-10 keV luminosity, L2-i0keV"^ 

1042 _ io44 ergs s-'AGNs). On average, the infrared SED of 
typical AGNs is best described as a broken power-law at < 40 |^m that falls steeply at >40 i_im 
(i.e., at far-infrared wavelengths). Despite this fall-off at long wavelengths, at least 3 of the 1 1 
AGNs in our sample have observed SEDs that are AGN-dominated even at 60 i.im, demon- 
strating the importance of accounting for possible AGN contribution even at far-infrared 
wavelengths. Our results also suggest that the average intrinsic AGN 6-100 ^.m SED gets 
bluer with increasing X-ray luminosity, a trend seen both within our sample and also when we 
compare against the intrinsic SEDs of more luminous quasars (i.e., L2-i0keV^ 10'*'* ergs s^*). 
We compare our intrinsic AGN SEDs with predictions from dusty torus models and find they 
are more closely matched by clumpy, rather than continuous, torus models. Next, we use our 
intrinsic AGN SEDs to define a set of correction factors to convert either monochromatic in- 
frared or X-ray luminosities into total intrinsic AGN infrared (i.e., 8-1000 [im) luminosities. 
Finally, we outline a procedure that uses our newly defined intrinsic AGN infrared SEDs, in 
conjunction with a selection of host-galaxy templates, to fit the infrared photometry of com- 
posite galaxies and measure the AGN contribution to their total infrared output. We verify the 
accuracy of our SED fitting procedure by comparing our results to two independent measures 
of AGN contribution: (1)12 |a,m luminosities obtained from high-spatial resolution observa- 
tions of nearby galaxies and (2) the equivalent width of the 1 1 .25 [im PAH feature. Our SED 
fitting procedure opens up the possibility of measuring the intrinsic AGN luminosities of large 
numbers of galaxies with well-sampled infrared data (e.g., IRAS, ISO, Spitzer and Herschel). 
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1 INTRODUCTION 

The spectral energy distribution (hereafter, SED) of a continuum 
source is the description of its energy output as a function of pho- 
ton frequency or wavelength. As such, it is one of the most impor- 
tant measurables in astronomy, providing information on both the 
physical nature of a continuum source (e.g., stars, galaxies, active 
galactic nuclei, heated dust), its influence on the surrounding matter 
(i.e., heating, ionisation state) and, when integrated over all wave- 
lengths, its bolometric luminosity (i.e., power output). However, 



* The intrinsic AGN infrared SEDs, host-galaxy templates and a pro- 
cedure used to combine them to fit infrared photometry is available at 
http://sites.google.com/site/decompir 
t E-mail: james.mullaney@cea.fr 



deriving an SED is a challenging task, requiring multiple observa- 
tions across the whole of the electromagnetic spectrum and, in the 
case of active galactic nuclei (hereafter, AGNs), one that is further 
complicated by contamination from the ever-present host galaxy. 

The level of host-galaxy contamination to the observed AGN 
SED is typically highest at infrared wavelengths (i.e., 8-1000 \im) 
whe re the strongly rising host-galaxy SED typically domin ates 
(e.g. lElvis et al.lll994l : lRichards et al.ll200dlNetzer et al.l2007h . In- 
deed, this contamination is so severe that the intrinsic AGN SED at 
these wavelengths remains largely unconstrained by observations 
(except for rare quasar-luminosity AGNs; see later). The signifi- 
cance of the uncertainties surrounding the intrinsic AGN infrared 
SED is clear when we consider that the infrared portion of the 
average observed (i.e., AGN + host) AGN SE D represents >5 
per cent of the total radio - X-ray output (e.g., IeMs et anil994h . 
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Physical models of the infrared-e mitting dust surrounding the AGN 
can provid e som e insights (e.g.. ISiebenmorgen &Kruegel|[T993 : 
iFritz et al.l I2OO6I: ISchartmann et al.l 12008). although they tend to 
predict a very broad range of infrared AGN SEDs as their input 
parameters are typically poorly constrained (see our Sj53J; a situa- 
tion that will be helped with a better understanding of the true range 
of intrinsic AGN infrared SEDs. 

Well-defined intrinsic AGN infrared SEDs may also prove a 
useful tool in the measurements of AGN and star-formation activ- 
ity that are crucial to our understanding of the inte ractions between 
these two processes that are believed to exist (e.g., Magorrian e t al 
1998; Ferrarese & Merritt 2000; Gebhardt et al. 2000; Croton et_al 
20061 ; iBower et alj 120061 ; iHopkins et al.l |2006|. 120071) . It has al- 



ready been demonstrated that the intrinsic mid-infrared emission 
of AGNs is clo s ely linked to the total AGN luminosity (e.g., 
iHorst et alj|2008l ; iGandhi et al.ll2009h , while the far-infrared emis- 
sion of a pure star-forming galaxy provides a proxy measure of 
its star formation rate (e.g.. lKennicutl|[l998h . Separating the total 
infrared SED of a composite (i.e., AGN + star-forming) galaxy 
into AGN and host components would therefore provide a mea- 
sure of the levels of AGN and star-formation activity that is largely 
unaffected by obscuration or absorption; similar, in principle, to 
the infrared spectral decompositions c arrie d out by, for exa mple, 
iLaurent et all ( l2000l)lTran et al.l ( 1200 lb and lLutz et'S] (|2004 |). The 
benefit of using the whole infrared SED is that this technique is 
not restricted to the small fraction of galaxies for which infrared 
spectra are available. Furthermore, such techniques are likely to 
become increasingly popular with the availability of infrared pho- 
tometry measurements from the Spitzer and Herschel telescopes. 
However, most studies that have used SED decomposition ap- 
proaches have had to rely on either observed AGN SEDs that 
likely include considerable amounts of host-galaxy contamination 
at infrared wavelengths or predictions from ra diative transfer mod- 
els which are often poorly constrained (e.g. Polletta et al. '2007"; 
Fioreetal. 2008, 2009; Pozzi et al. 2010; Hatziminaoglou et al. 
2010) . Of course, measurements derived from decomposing the in- 



frared SED should be regarded as complementary to other mea- 
sures of AGN activity that rely on other observables (e.g.. X-ray lu- 
minositVj^ optical and infrared spectroscopy, mid-infrared photome- 
try etc.;lKau fj mann_ el^aL 2003; Heckman et al. 2004; Hickox et al. 
l200ll2009l ; lGoulding & Alexanderil2 009). 

In their I2OO7I paper, iNetzer et a l. derived the average intrin- 
sic infrared SED of a sample of luminous PG quasars by sub- 
tracting the host-galaxy component from the average observed 
SED, showing that their average intrinsic AGN SED falls rapidly 
longward of ~ 20 \im. However, such luminous quasars are rare 
among the AGN population and it is not clear whether the in- 
trinsic quasar SED can be applied to the much larger, less lu- 
minous (i.e., L2-l0keV < 10^^ ergs s^') population of AGNs 
which represents the majorit y of the integrated accretion output of 
AGNs in the Universe (e.g.. lUeda et al. 



jgrated accretion output or 
l l2003llBarger et al.ll2005l ; 



iHasinger et al.' '2005*; Aird et al. '2010*). In this paper we extend 
the work of Netzer et al. (2007) to include more typical, lower lu- 
minosity AGNs and demonstrate that the intrinsic AGN infrared 
SED appears to be linked to AGN luminosity. To do this, we first 
identify a sample of AGNs that show minimal amounts of host- 
galaxy contamination at MIR wavelengths, then carefully subtract 
any host-galaxy component that can still dominate at far-infrared 
wavelengths (see ij2] Sj3]l. We compare the derived intrinsic AGN 
SEDs against other, more commonly assumed, AGN SEDs includ- 
ing those derived from quasars and radiative transfer models in ^ 
In ^we introduce a procedure that uses these intrinsic AGN SEDs, 
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Figure 1. The distribution of intrinsic 2-10 keV luminosities (i.e., L2-I0kev) 
and hydrogen column densities (i.e., Nn) of the Swift-BPCt sample of AGNs 
which we use to de fine the intrinsic AGN infrared SED (circles; taken from 
IWinter et al.ll200^ . In the background (small squares) we also show the 
^2-lOke V~Mt distributio n of galaxies in the Chandra Deep Field-South (i.e., 
CDF-S; Tozzi et al. 2006). The intrinsic properties of these local (i.e., Swift- 
BAT) and high-redshift (i.e., CDF-S) samples of galaxies ai'e well matched. 
We highlight those 25 and 1 1 AGNs (see key) that we use to define the 
range of intrinsic AGN mid-infrared and mid to far-infrared SEDs. These 
subsamples cover a similai' range of i2-lOkeV^Mt parameter space spanned 
by the Swift-BPiT and CDF-S samples. See online manuscript for a colour 
version of this plot. 



together with a sample host-galaxy templates, to measure the AGN 
and host-galaxy contributions to the infrared output of a sample of 
local composite galaxies and demonstrate the obtained values agree 
with results from other, independent measures. In ^Zjwe outline two 
obvious applications of these analyses: (1) defining correction fac- 
tors to convert 12 (im and X-ray luminosities to total intrinsic AGN 
infrared luminosities and (2) measuring the intrinsic AGN power 
using infrared photometry alone. We summarise our findings in ^ 
Throughout this work we adopt Hq = 11 km s^' Mpc^', 
flM = 0.27, and Slj^ = 0.73. 



2 SAMPLE DESCRIPTION 

The first step we take in defining the intrinsic infrared SEDs 
of typical AGNs is to identify a sample of moderate luminosity 
(i.e., i-a-lOkeV^ 10^^ ^ 10^^ ergs s^') AGNs that suffer minimal 
amounts of host galaxy contamination at infrared wavelengths. By 
doing this, we only focus on those SEDs that require as little ma- 
nipulation as possible to extract the intrinsic AGN infrared SED. In 
this section, we describe how combine information derived from 
hard X-ray observations and infrared spectra to produce such a 
sample. We then explain how we use IRAS photometry to extrapo- 
late the SEDs of a subsample of these AGNs to FIR wavelengths. 

The Burst Alert Telescope (hereafter, BAT) on-board the Swift 
observatory is currently undertaking a survey of the sky at hard 
X-ray energies (i.e., 14-195 keV), where the effects from ab- 
sorption (at least to A'h< lO^'* cm^^) are negligible. As a re- 
sult, the second data release of the 5wi/f-BAT survey, published 



in iTueller et al. provides a homogeneous sample of 154 

local (i.e., z < 0. 1) X-ray AGNs. A key motivation for using the 
AGNs selected from the Swift-BAT sample to define the intrinsic 
AGN infrared SED is that their X-ray properties (i.e., 2-10 keV X- 
ray luminosities and absorbing columns, hereafter, L2_iokeV ™d 
A'h) cover largely the same ranges as those AGNs detected in 
deep Chandra surveys (e.g., CDF-N and CDF-S; see fig. [B- As 
such, this sample provides the ideal local analogue to the high red- 
shift AGNs detected in these fields as well as typical local AGNs 
(i.e., L7.-iofav~ 1 0"^^ - 10'^'* ergs s"' and Nii< lO^"* cm^^. 
IWinter et al.i2009[) . Furthermore in cases where infrared photom- 
etry is available, the Swift-BAT AGNs cover the same range of 
infrared luminosities (integrated over 8-1000 (xm; hereafter Lir; 
LiR= 6 X 10^ — 7 X 10" L,7,) as those star-forming galaxies that 
make up the majority of the local galaxy population detected by 
IRAS (i.e., LiR~ 10'° - lO'^ Lq) and z < 2 galaxies detected in 
deep SpitzerlHerschel fields (e.g. GOODS; P.I.s: M. Dickinson, D. 
Elbaz). 

In order to identify a sample of galaxies whose infrared out- 
put is dominated by a moderate luminosity AGN, we cross-match 
the 5wi//-BAT sample with the archive of low resolution spectra 
obtained by the infrared spectrograph (hereafter, IRS) on-board 
Spit zer. Of the 104 Sw ift-BAT AGNs with measured L2_iokeV 
(see I Winter et al.ll2009l) . 36 have publicly available archival low 
resolution IRS spectra covering (observed frame) '--^6 — 35 yum. 
These spectra were reduced a nd extracted following the proce- 
dures outlined in lMuUanev et al. (2010), Goulding et a l. ( 201(1) and 
G oulding (2010; PhD thesis). Fo llowi ng the diagnostic s prese nted 
in iGoulding & Alexander! ( |2009|) and iTommasin et al.l ( |2010|) . we 
quantify the contribution of the host galaxy to these MIR spectra 
using the equivalent widths of the 11.25 [im PAH feature (here- 
after, EW PAHj^ji 25). We decide to use the 11.25 [im feature 
to measure the host galaxy contribution rather than the typically 
stronger one at 7.7 \im as this latter feature is often blended with 
another PAH emission line at 8.6 [im. We measure EW PA H^j j 1 75 
using the spectral fitting code PAHFIT jSmith et al]|2007h and as- 
sume that strongly AGN dominated systems have EW PAH^j^n 25 < 
0.03 [im, corresp onding to a < 10 per c ent host-galaxy contribu- 
tion at 19 |i.m (see' Tonimasin et alfcOlOh Fl Of the 36 X-ray AGNs 
that have archival low resolution IRS spectra, we identify 25 whose 
EW PAH^^i 1 25 satisfy this criterion (see table[T|for a list of these 25 
AGN-dominated galaxies; also fig.|2j. Approximately half of these 
AGNs (i.e., 12/25) are optically classified as either Type 1, 1.2 or 
1.5 while the rest (i.e., 13/25) are classified at either Type 1.8, 1.9 
or 2 (optical classifications taken from Tueller et al. 2008; see table 
[T]for these classifications and other information regarding these 25 
AGN dominated galaxies). We visually inspect the infrared spectra 
of each of these 25 AGNs to confirm that they are, indeed, domi- 
nated by a featureless AGN continuum. 

To provide constraints on the intrinsic AGN SED at FIR wave- 
lengths, we extrapolate beyond the low resolution IRS spectra using 
60 vim and 100 (im photometry available from the IRAS archives 
(accessed through the NASA extragalactic database; NED). To en- 
sure reliable extrapolation to FIR wavelengths, we only consider 
photometry greater than 0.2 ly and 1.0 ly at 60 [im and 100 yim, 
respectively (the approximate limiting sensitivity of the IRAS faint 
source catalogue at these wavelengths). Of the 25 AGN-dominated 
galaxies identified above, 20 satisfy these criteria. However, we 
note that aperture effects must be taken into consideration when 



collating photometric and spectral data obtained by the IRAS and 
Spitzer telescopes. The apertures used to measure the IRAS pho- 
tometry reported in the point source and faint source catalogues 
were between 2' and 9' in diameter, compared to the 3.6" and 10.5" 
(shortest dimension) apertures of the short-low and long-low slits of 
the Spitzer-TRS instrument, respectively. Therefore, there is consid- 
erable scope for the IRAS photometry to contain a great deal more 
host galaxy flux than the small aperture IRS spectra. To mitigate 
such aperture effects, we only extrapolate to the 60 [im and 100 [im 
photometry values when the flux density in the IRS spectrum, in- 
tegrated over the 12 and 25 [im IRAS band passes, matches (to 
within 30 per cent, including eiTors) the IRAS photometry at these 
wavelengths. Our choice of 30 per cent accuracy between the IRAS 
and IRS fluxes is based on the approximate systematic uncertainty 
of the IRAS point source and faint source c atalogues when com - 
pared to the Revised Bright Galaxy Sample jSanders et al.ll200^ . 
We identify 11 AGN-dominated galaxies that satisfy these crite- 
ria (indicated by a "Y" in column 16 of table [T), which we use 
to define the intrinsic AGN SED. The 6-100 nm SEDs of these 
AGN-dominated galaxies, including their 60 \im and 100 [im pho- 
tometries, are presented in appendix [Al 



PAHFIT is available from |http://tir.as tro.utoledo.edu/jdsmith/research/pahfit.php | 
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Figure 2. The 25 AGNs from the Swift-JiAI sample that are strongly AGN dominated at mid-infrared wavelengths, separated in terms of their optical 
classification (here. Types 1.2 and 1.5 are classed as Type 1 AGN, Types 1.8 and 1.9 are classed as Type 2 AGNs). We have indicated the positions of the 
most common, prominent emission lines, PAH features and silicate emission/absorption features. The criterion we use to determine whether the mid-infrared 
SED is AGN-dominated i s that the equivalent width of the 1 1.25 |j.m PAH feature is less than 0.03 \un, which is equivalent to < 10 per cent host-galaxy 
contribution at 19 \xm (e.g. lTommasin et alj201 0). Within this subsample we see a wide range of mid-infrared SED shapes. However, in the majority of cases 
(i.e., at least 20/25; asterisked) the underlying AGN continuum can be described as an absorbed broken power law with a higher (i.e., more positive) spectral 
index longward of a break at, roughly, 19 |xm. In general. Type 2 AGNs have steeper SEDs (i.e., high spectral indices) at A, < 19 |xm although, in this respect, 
there is considerable overlap between the two classes. Type 2 AGNs are more likely to show evidence of strong silicate absorption at ~10 |J.m, while a larger 
proportion of Type 1 AGNs have silicate in emission at ~10 [im and ~18 |a.m. The galaxy names have been truncated to enable them to fit on the plot. See 
table[T]for the full galaxy names. 



Table 1. Properties of the 25 AGNs in the Sw(//-BAT sample that are strongly AGN dominated at 6-35 (J.m. 



Name 


Type 


RA 


Dec. 


Dist. 


Sii 


Sis 


^60 


5ioo 


cIRS 


eIRS 


vLv(60|j.m) 


^-IR 


i'2-lOkcV 


EWPAHxu.25 


Selected 










(Mpc) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


log(L,i) 


log(i0) 


log(ergs s ) 


(_^^m) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


Mrk348 


2 


12.1964 


31.9570 


63.4 


0.308 ± 0.031 


0.835 ± 0.025 


1.29 ± 0.12 


1.55 ± 0.20 


0.1 


0.4 


9.91 


10.5 


42.6 


0.002 




NGC985 


1 


38.6574 


-8.7876 


185.3 


0.207 ± 0.031 


0.523 ± 0.017 


1.38 ± 0.11 


1.89 ± 0.15 


0.1 


0.3 


10.87 


11.3 


44.3 


0.016 




NGC1275 


2 


49.9507 


41.51 17 


74.2 


1.060 ± 0.014 


3.440 ± 0.027 


6.990 ± 0.042 


7.20 ± 0.47 


0.7 


2.8 


10.78 


11.2 


43.9 






ESO548-G081 


1 


55.5155 


-21.2444 


61.0 


0.249 ± 0.057 


0.097 ± 0.075 


0.601 ± 0.048 


1.49 ± 0.18 


0.2 


0.3 


9.54 


10.2 


42.8 


0.015 




3C120 


1 


68.2962 


5.3543 


141.2 


0.289 ± 0.043 


0.624 ± 0.094 


1.38 ± 0.21 


1.94 ± 0.29 


0.3 


0.8 


10.63 


11.1 


44.0 


0.008 


Y 


NGC2110 


2 


88.0474 


-7.4562 


32.6 


0.349 ± 0.028 


0.840 ± 0.021 


4.13 ± 0.21 


5.7 ± 1.4 


0.3 


0.8 


9.83 


10.2 


42.5 


0.023 


Y 


EXO055620-3820.2 


1 


89.5083 


-38.3346 


144.9 


0.529 ± 0.032 


0.685 ± 0.034 


0.322 ± 0.035 


0.56 ± 0.12 


0.3 


0.6 


10.02 


11.2 


43.7 






Mrk3 


2 


93.9015 


71.0375 


56.9 


0.713 ± 0.050 


2.900 ± 0.036 


3.77 ± 0.15 


3.36 ± 0.44 


0.6 


2.3 


10.28 


10.8 


42.9 




Y 


MCG-0 1-24-0 12 


2 


140.1927 


-8.0561 


83.1 


0.080 ± 0.025 


0.3810 ± 0.0088 


0.654 ± 0.059 


1.10 ± 0.18 


0.1 


0.4 


9.85 


10.4 


43.0 


0.006 


Y 


Mrk417 


2 


162.3789 


22.9644 


140.1 


0.132 ± 0.033 


0.227 ± 0.054 


0.164 ± 0.043 


0.71 ± 0.15 


0.1 


0.3 


9.70 


10.7 


43.3 


0.004 




NGC3783 


1 


174.7572 


-37.7386 


38.5' 


0.840 ± 0.059 


2.490 ± 0.050 


3.26 ± 0.19 


4.90 ± 0.54 


0.6 


1.9 


9.88 


10.5 


43.1 


0.002 




NGC4507 


2 


188.9026 


-39.9093 


49.6 


0.517 ± 0.078 


1.59 ± 0.24 


4.69 ± 0.70 


6.28 ± 0.94 


0.6 


1.5 


10.25 


10.7 


43.0 


0.013 


Y 


ESO506-G027 


2 


189.7275 


-27.3078 


106.3 


0.148 ± 0.030 


0.268 ± 0.022 


0.498 ± 0.050 


0.83 ± 0.19 


0.2 


0.4 


9.94 


10.5 


43.2 






MCG-03-34-064 


1.8 


200.6019 


-16.7286 


69.8 


0.940 ± 0.040 


2.970 ± 0.045 


6.200 ± 0.040 


6.20 ± 0.14 


0.7 


2.5 


10.67 


11.1 


42.8 


0.005 


Y 


MCG-06-30-015 


1.2 


203.9741 


-34.2956 


32.5 


0.380 ± 0.034 


0.809 ± 0.027 


1.090 ± 0.076 


1.10 ± 0.22 


0.2 


0.5 


9.25 


9.9 


42.6 


0.004 




IC4329A 


1.2 


207.3304 


-30.3096 


67.7 


1.080 ± 0.054 


2.210 ± 0.054 


2.03 ± 0.10 


1.66 ± 0.22 


0.9 


2.0 


10.16 


10.9 


43.6 




Y 


NGC5506 


1.9 


213.3119 


-3.2075 


28.7' 


1.290 ± 0.027 


4.170 ± 0.056 


8.420 ± 0.060 


8.87 ±0.11 


1.0 


3.3 


10.03 


10.5 


42.8 


0.011 


Y 


NGC5548 


1.5 


214.4981 


25.1368 


72.5 


0.401 ± 0.040 


0.769 ± 0.032 


1.070 ±0.086 


1.61 ±0.16 


0.2 


0.6 


9.94 


10.6 


43.5 


0.017 




Mrk290 


1 


233.9682 


57.9026 


126.1 


<0.1 


<0.1 


0.171 ±0.029 


<0.6 


0.1 


0.2 


9.63 


<10.5 


43.2 






ESO103-035 


2 


279.5847 


-65.4276 


55.9 


0.612 ± 0.043 


2.360 ±0.031 


2.31 ±0.12 


1.05 ±0.26 


0.5 


2.1 


10.05 


10.7 


42.7 




Y 


3C390.3 


1 


280.5375 


79.7714 


244.3 


0.128 ±0.018 


0.2870 ± 0.0089 


0.204 ± 0.033 


<0.6 


0.1 


0.3 


10.28 


<11.2 


44.1 


0.002 




Mrk509 


1.2 


311.0406 


-10.7235 


147.3 


0.316 ±0.028 


0.702 ± 0.022 


1.360 ±0.068 


1.52 ±0.23 


0.3 


0.8 


10.66 


11.2 


44.0 


0.023 


Y 


IC5063 


2 


313.0097 


-57.0688 


47.7 


1.110 ±0.023 


3.940 ± 0.030 


5.870 ± 0.038 


4.25 ± 0.21 


0.9 


3.1 


10.32 


10.8 


42.7 


0.000 


Y 


NGC7213 


1.5 


332.3177 


-47.1667 


22.0' 


0.606 ± 0.049 


0.742 ± 0.036 


2.67 ±0.16 


8.18 ±0.41 


0.2 


0.5 


9.30 


9.9 


42.3 


0.019 




NGC7314 


1.9 


338.9426 


-26.0502 


19.0' 


0.268 ± 0.035 


0.579 ± 0.048 


3.74 ± 0.22 


14.2 ± 1.3 


0.1 


0.5 


9.32 


9.8 


42.2 


0.015 





Notes: (1) Common name, (2) Optical Class from lTueller et al.l j2008h (3) and (4) Galaxy co-or dinates (J2000) from lTueller et al.l j2008l) . (5) Distance in Mpc, 'denotes that a redshift-independent distance measure 
taken from NED is used, otherwise calculated from the redshift reported in iTueller et al J i2008l) (6) to (9) IRAS photo metry measurement s, taken from NED (10) and (11) 12 |j.m and 25 [im flux densities derived 
from the low resolution IRS spectra, (12) Monochromati c luminosity at 60 \ i m in u nits of log(L0) (for comparison with lNetzer et al .l2007h . We use L0= 3.83 X 10^^ ergs s '. (13) Total infrared luminosity in units 
of log(L;,), calculated using the equations in Table 1 of lSanders & Mrrabell jl996l) . (14) Rest-frame 2-10 keV luminosity, corrected for absorption (taken from Iwinter et al.ll200 9^ (15) Equivalent width of the the 
1 1.25 |J.m PAH (i.e., EW PAH^^n 25) feature in |j.m. (16) Selected for extrapolation to 100 [im based on whether the IRAS photometry at 12 and 25 |J.m matches that derived from the IRS spectra. 



3 DEFINING A SET OF HOST-GALAXY TEMPLATES 

The sample of AGNs identified in the previous section show 
minimal amounts of host-galaxy contamination at MIR wave- 
lengths. However, because a typical host-galaxy SED rises strongly 
at longer infrared wavelengths, an SED that is strongly AGN- 
dominated at MIR wavelengths may not necessarily b e AGN- 
dominated at FIR wavelengths. Indeed, iNetzer et ail ( l2007h showed 
that at 60 (xm the average quasar SED is 80 to 90 per cent host- 
galaxy dominated, despite having only weak PAH features at MIR 
wavelengths. In this section we describe the construction of a 
set of host-galaxy templates that represent the full diversity of 
host-galaxy SEDs expected for typical AGNs (i.e., the full range 
of IRAS colour-colour space for LiR= 10'° - 10'2 Lq galaxies). 
These templates are used to account for the host-galaxy contribu- 
tion to our AGN-dominated sample in jj?) 

In fig. [3] (left panel) we plot the [100 (j.m/60 \un\ vs. 
[25 jxm/12 i^m] colours of all Lir= I0'° — lO'^ galaxies in 
the Revised Bright Galaxy Survey. From this figure it is clear that 
the majority of the colou r-colour spac e span ned by these galax- 
ies is well sampled by the lSrandl et alj ( |2006|) sample of starburst 
galaxies. However, we note that there is a region of this colour- 
colour space that is not sampled by these starbursts. Galaxies in this 
un-sampled region of the colour-colour space have SEDs that rise 
steeply between 60 ]mi and 100 (xm, but are relatively flat between 
12 \im and 25 \im. As such, they likely represent a population of 
cold gala xies. To construct ou r host-galaxy templates we use a se- 
lection o f lBrandletal.1 ( l2006h starbursts and a selection of galaxies 
that lie in this "cold" region of the IRAS colour-colour space. This 
ensures that we cover the full range of SEDs displayed by Lir— 
lO'" — lO'^ Lp: galaxies. We note that the both the colour-colour 
and [60 ^m/25 jim]-LiR parameter spaces are particularly well- 
sampled in the infrared luminosity range spanned by the AGNs se- 
lected for decomposition (i.e., \og(Li^/L(r,)= 10.2— 1 1.2 Lp,; fig.[3] 
right panel; see ^4.21 . 

We extracted an d reduced the low resolution IRS spectra of all 
the lBrandl et alj ( [2006) starbursts and four additional "cold" galax- 
ies (namely, NGC 1667, NGC 5734, NGC 6286 and NGC 7590), 
following the same procedures as those referred to in ij2] All of the 
iBrandletalJ JIoOe) starbursts were observed in IRS staring mode, 
while the four "cold" galaxies were observed in mapping mode. We 
note that these differences in the observing modes will have little or 
no effect on our host-galaxy templates as our selection process min- 
imises aperture effects (see later). We extrapolated these MIR SEDs 
to 100 ^.m using IRAS photometry as described in ij2] To mitigate 
aperture effects we only use those SEDs in which the IRS spectra 
agree with the 12 [i m and 25 [im IRAS p hotometry to within 30 per 
cent. Out of the I6 lBrandl et a starbursts identified as ei- 

ther pure starbursts or starburst +LINER, we identify 10 that satisfy 
these criteria: Mrk 52, Mrk 520, NGC 660, NGC 1222, NGC 2623, 
NGC 3256, NGC 4194, NGC 4818, NGC 7252, NGC 7714. All 
four of the "cold" galaxies also satisfy these criteria. To produce 
our host-galaxy templates we group these SEDs in terms of their 
overall shape and the relative strength of their PAH features (see fig. 
|4ll. The full range of host-galaxy SED shapes is well characterised 
by five groups of SEDs, referred to as "SBl" through "SB5", with 
all four "cold" SEDs contained within a single group (i.e., "SB I"; 
see table [2] and fig. @. Five host-galaxy templates are obtained by 
noiTnalising each SED within each group at 90 [im and calculating 
their mean average. Since the SEDs within each group are so sim- 
ilar, the wavelength at which they are normalised has little effect 
on these host-galaxy templates. Beyond 100 [im we extrapolate the 



average SEDs as a modified black body (i.e., Fv = E^^v^, where 
F^^ is the blackbody specific flux and v is photon frequency; we 
adopt (3 = 1.5). We verify that this is a reasonable assumption using 
AKARI 140 |i.m and 160 vim photometry data available for 13 of 
the 14 host-galaxy SEDs in our sample^ These templates are used 
to remove any host-galaxy contribution to the infrared SED of the 
AGN-dominated galaxy sample. We publish the host-galaxy tem- 
plates in columns 5-9 of table [3] (which is available in its entirety 
online at ,http://sites. google.com/site/decompirj . 



4 CHARACTERISING THE INTRINSIC INFRARED AGN 
SED 

First we use the IRS spectra of our AGN sample (defined in §2.1) 
to explore the range of AGN-dominated infrared SEDs at 6-35 vim. 
We then extract the intrinsic 6-100 [im SEDs of a subsample of 
these AGNs by carefully decomposing the observed SEDs into their 
host-galaxy (defined in ij3j and AGN components. 

4.1 AGN dominated SEDs at 6 vim to 35 vim 

The IRS spectra presented in fig. |2] provide a clear picture of the 
range of MIR SEDs produced by strongly AGN-dominated sys- 
tems. As explained in ij2l the weak PAH features in this sample 
suggest that at least 90 per cent of the continuum emission at 19 vim 
is produced by the AGN. Considering that the typical host-galaxy 
continuum falls toward shorter infrared wavelengths (see fig. |4]l, 
it is reasonable to assume that the AGN dominates the contin- 
uum emission at ~6 vim - 20 vim. At these wavelengths, all 25 
AGN-dominated SEDs show clear evidence of a continuous, un- 
derlying power-law continuum that is thought to be produced by 
multiple dust compon ents spanning a range of temperatures (e.g., 
[Buchanan et"aill2006h . 

The spectral indices of the underlying AGN power-law con- 
tinua span the range 0.7 ^ ai ^ 2.7 (mean: a; = 1.6), with a ten- 
dency for Type 2 AGNs to have higher (more positive) spectral in- 
dices compared to Type 1 AGNs. |j However, there is considerable 
overlap between the types (i.e.. Type Is: 0.7 ^ ai ^ 1.7; Type 2s: 
0.8 ^ CXi ^ 2.7), and the difference in the mean SEDs is not sta- 
tistically significant. This large range of spe ctral indices is consis- 
tent with those found in previous studies (e.g. Buchanan et al.'2006|; 
IWu et alj2069h . In four cases (i.e., one Type 1 : EXO055620-3820.2 
and three type 2s: ESO506-G027, NGC 5506, ESOI03-035), we 
see evidence of strong silicate absoiption at 9.7 vim. In a further 
eight cases (i.e., five type Is: ESO548-G081, 3C120, NGC 3783, 
Mrk 290, NGC 7213 and three type 2s: NGC 1275, NGC 2110, 
NGC 4507) there is strong evidence of silicate emission at ^ 10 vim 
which is always accompanied by a nother silicate emission feature 
at ^ 18 vim (see lSturm et al .I2OO5 I for a dedicated study of the sili- 
cate emission features seen in the MIR spectra of AGNs). 

In at least 20 of the 25 AGN-dominated MIR spectra there is 
a definite break in the power-law continuum at 15 — 20 vim; mean 
break position of ~ 19 vim. In previous studies this has been at- 
tributed to a dominating wa rm (i.e., 170 K ) dust compon ent that 
is heated by the AGN (e.g.. lWeedman et"al.ii2005i : .Buchanan et al.l 

^ AKARI data were retrieved from the NASA/IPAC Infraied Science 
Archive at http://rrsa.ipac.caltech.edu/ 

' Here, we use /\^"'^ « A.". To convert to v/\^™ use (3 = a — 1, where 
Vf AGN „ ^f, 
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Figure 3. IRAS colour-colour (left panel) and colour-luminosity (right panel) plots of all sources in the Revised Bright Galaxy Sun'ey jSanders et alj|2003[) 
with infrared luminosities covering the roughly the same range as the Swift- BAI sample from which our AGN-dominated galaxies were selected (i.e., Lir= 
lO'" — lO'^ ergs s"'). We use a subsample of 10 starburst SEDs from Brandl et al. (2006) plus four other "cold" galaxies to define our host galaxy templates 
(circled stars). These extra four q uiescent-galax y SEDs are needed to sample the full range of host-galaxy SEDs covered by Lir= lO'" — lO'^ ergs s"' galaxies. 
The lack of "cold" galaxies in the lBrandl et alj (2006) sample is shown most clearly in the colour-colour plot (left hand panel). The horizontal black line in the 
right hand panel indicates the range of Lir of the AGNs selected for decomposition (see i]4.2t . The objects with [60 |j.m/25 |xm] < 4 are likely AGN dominated 
and were therefore not considered as suitable host-galaxy templates. 



In all cases where we see a break in the continuum power 
law, the AGN SED longward of X^a has a lower spectral index 
(i.e., po wer-law index, 0^ go ^ 1.5, mean= 0.7 ; again, consis- 
tent with lBuchanan et al.ll2006f and lWu et alj|2009l) than at shorter 
wavelengths. We again find that there is significant overlap in the 
range of spectral indices at X> 19 \im between Type 1 and Type 
2 AGNs. We note that the power-law indices shortward and long- 
ward of ^Brk ^6 uncorrelated. The number of AGNs in our sample 
showing a break around 19 [im could be as high as 24, since it 
can be disguised by the presence of the silicate emission feature at 
18 [im. The only AGN in our sample where the 6-35 [im SED can 
be unambiguously described as a single power law is NGC5506; 
the spectral index of NGC5506 over 6-35 \Lm is a = 1.2 ±0.2, 
which is consistent with the ranges of both ai and a2- 

Because the typical host-galaxy continuum emission increases 
strongly toward longer wavelengths, it is difficult to ascertain 
whether the power-law emission at longer wavelengths (i.e., A. > 
25 \im) arises from the AGN alone, or if the host galaxy also makes 
a significant contribution. Indeed, in four cases (i.e., NGC7213, 
NGC4507, NGC2110, 3C120), the SED longward of ^ 25 ^m 
shows evidence of a turn-up in vFy which is consistent with be- 
ing due to an underlying host-galaxy component. The host-galaxy 
component must be accounted for at longer wavelengths in order 
to define the intrinsic AGN infrared SED, which we address in the 
next subsection. 



4.2 Intrinsic AGN infrared SED at 6 ^m to 100 

We now extend the intrinsic AGN infrared SEDs to longer wave- 
lengths (i.e., X>25 [im) using IRAS 60 [im and 100 yum photome- 



try for guidance. Here, we only consider those 1 1 AGN-dominated 
SEDs that we have shown suffer from minimal aperture effects (see 
The approach we take is to decompose the observed SEDs into 
their host-galaxy and intrinsic AGN components through simulta- 
neously fitting the IRS spectra and IRAS photometry. An alterna- 
tive approach would be to simply subtract a normalised host-galaxy 
component from the observed SEDs to leave the emission intrinsic 
to the AGN. However, calculating the appropriate normalisation for 
the host-galaxy component is a non-trivial matter that is compli- 
cated by the fact that the host-galaxy SED may be modified by 
dust-absorption (i.e., the observed SED is not a simple linear sum 
of a host-galaxy and an AGN component; separate absorption terms 
must also be considered). Later, and only for illustrative purposes, 
we subtract the host-galaxy SEDs normalised via our SED fitting 
routine from the sample of 1 1 observed AGN SEDs (see fig.H). 

One downside of the approach that we adopt is that it requires 
us to make a priori assumptions about the form of the intrinsic 
AGN infrared SED used in our fits. Our analyses of the AGN- 
dominated MIR spectra in the previous subsection provides us with 
a well defined set of parameters that describe the variety of intrinsic 
AGN SEDs at 6 to ~25 ^m (i.e., 0.7 < ai ^ 2.7, < < 1 .5 and 
15 nm^ Xg^^ait^ ^20 \im). At longer wavelengths, where there are 
fewer observational constraints, we can only estimate the form of 
the intrinsic SED based on reasonable physical assumptions then 
test whether models incorporating these SEDs reproduce observa- 
tions. With this in mind, we assume that the intrinsic SED falls as 
a modified black body beyond a given wavelength, A-bb- In our fits, 
A-BB is allowed to take any value between 20 [im and 100 [im. Our 
choice of a modified black body spectrum is loosely based on the 
results of radiative transfer models of the dust surrounding AGNs, 



Table 2. Galaxy SEDs that used to construct our host galaxy templates. 



10' 



Host Galaxy Template 


Galaxy 


D (Mpc) 


log(LiR/L0) 


(1) 


(2) 


(3) 


(4) 


SBl 




60.51 


10 96 




NGC 5734 


59.28 


11.06 




NGC 6286 


79.78 


11.32 




NGC 7590 


21.58 


10.16 


SB2 


NGC 7252 


64.67° 


10.77'' 


SB3 


Mrk52 


33.50* 


10.25^ 




NGC 4818 


9.37 


9.75 




NGC 7714 


38.16 


10.72 


SB4 


NGC 1222 


32.26 


10.60 




NGC 3256 


35.35 


11.56 




NGC 4194 


40.33 


11.06 


SB5 


NGC 520 


30.22 


10.91 




NGC 660 


12.33 


10.49 




NGC 2623 


77.43 


11.54 



Notes: (1) Host-galaxy templat e, (2) Galaxy name, (3) Distance in Mpc 
taken from ISanders et al.l j20O3l) . "calculated using the redshift reported 
in NED using the same cosmology as Sanders et al. (2003), ''redshift- 
independen t distance reported in NED (4) 8-1000 \ua infrared luminosity, 
taken from lSanders et all l l2003h . '^calculat ed using the distanc e in co lumn 
3 and the prescription outlines in table 1 of lSanders & Mirabell jl996l) . 

although we note that the precise form of this fall-off does not have 
a significant impact on the fits. When needed, we add silicate emis- 
sion features at 10 |i.m and 18 (xm to the intrinsic AGN component. 
These are modelled with broad (i.e., FWHM~3 iJ.m) gaussians. In 
our models, we account for any absorption to both the AGN and 
host-galaxy components using a lDrainel ( [2003h extinction curve. 

We proceed by simultaneously fitting the IRS spectra and the 
IRAS photometry of the 1 1 AGNs in our sample with a combina- 
tion of our host-galaxy templates (see ij3} and the estimated intrin- 
sic AGN SED described above. Any strong emission lines (exclud- 
ing PAH features) are masked to ensure that they do not adversely 
effect the fit. Each of the 1 1 infrared SEDs are fit five times, each 
time using a different host-galaxy template (defined in ^ see also 
fig. Illl. We use x^-minimisation to obtain the best fitting parame- 
ters for each choice of host-galaxy template. Plots showing the fits 
to the data (including the host-galaxy and intrinsic AGN compo- 
nents) and the resulting residuals are presented in appendix IaI 

Because the observed SEDs consist, in part, of very high 
signal-to-noise IRS spectra, none of the fits to the data are formally 
good (i.e., their reduced y} » 1). Therefore, we cannot use 
statistics to unambiguously determine whether any of the five dif- 
ferent AGN~host-galaxy template solutions provide a good charac- 
terisation of the data. Instead, we use the following criteria to deter- 
mine which of the five models provide suitable fits to the observed 
data: does it (a) pass within 2c of the 60 |J.m and 100 |i.m photom- 
etry and (b) reproduce the general shape of the IRS spectrum. For 
each of the 1 1 SEDs in our sample, we can identify at least one of 
the five model fits that satisfy both these criteria. When both criteria 
are met we select, by eye, those model fits that produce the small- 
est residuals at the wavelengths of the strongest PAH features, i.e., 
7.7 (j.m, 1 1.25 \im and 12.8 vim. We use this additional criterion to 
ensure that any contribution from the host-galaxy well accounted- 
for by the fitted host-galaxy component. Where there is no signif- 
icant difference in the PAH residuals of two or more such fits we 
identify each of them as being suitable (in appendix |A] we have 
highlighted all selected fits). In only two cases, namely NGC 5506 
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Figure 4. The fiv e average host-galaxy templates derived from the 
iBrandl et al] <2006l) starburst sample and t he four "cold" galax ies selected 
from the Revised Bright Galaxy Sample jSanders et al]|2003l : see fig. |3]. 
Each template is the mean-average of a group of SEDs that have similar 
overall shape and relative PAH strengths (see [|3]for a description of how 
we construct these host-galaxy templates). Each SED is made up of a low- 
resolution IRS spectrum that we extrapolate to far-infrared wavelengths us- 
ing IRAS photoinetry at 60 n-in and 100 |j.m. We also show the mean flux at 
65 n-m, 140 n-in and 160 n-in measured from AKARI data (when available) 
as a check that the modified black-body extrapolation beyond 100 \im is a 
reasonable approximation of the real SEDs. We use these average templates 
to model any host galaxy components in our sample of AGN-dominated 
galaxies to deiive the intiinsic AGN infrared SED. 



and ESOI03-035, do we see large residuals at the PAH wavelengths 
for all five fits. In both these cases, we assume that all of the fits are 
as good as each other and select all five for further consideration. 
In all selected cases (including NGC 5506 and ESOI03-035), the 
model fits lie within 10 per cent of the measured flux density for 
at least 85 per cent of the IRS data points, demonstrating the abil- 
ity of our SED fitting approach to reproduce the continuum shape 
and PAH strengths of the AGNs in our sample. As a final note to 
the fitting procedure we would like to point out that, while we are 
confident that selecting the fits that closely reproduce the PAH fea- 
tures is the most appropriate procedure to take, our main results do 



not change if, instead, we use the resuhs derived from using all our 
host-galaxy templates (which cover the full range of SED shapes 
of local, LiR= 10^" — lO'^ Lq, galaxies; see fig. [3j. The strongest 
effect that taking this alternative approach has on our results is to 
increase the scatter of intrinsic AGN SEDs at FIR wavelengths by 
a factor of ~ 2 — 3; as we shall see, this does not represent a large 
increase over the range of intrinsic AGN SED shapes produced by 
selecting only the suitable host-galaxy components identified using 
the selection criteria outlined above. 

For illustrative purposes, we show in fig.|5]the SEDs obtained 
by subtracting the suitable host galaxy components (i.e., those ex- 
tracted from the suitable model fits identified above and highlighted 
in appendix|Al( from the observed SEDs. As we have selected only 
those fits that reproduce the general shape of the observed SED, the 
SEDs presented in fig.|5]are equivalent to the intrinsic AGN compo- 
nents produced during the fitting procedure and shown in appendix 
lAl For guidance, we have highlighted the SED that is produced by 
the fit with the smallest overall residuals between 6 (im and 35 [J-m 
(i.e., the "best-of-the-best"), although we stress that all plotted in- 
trinsic SEDs are viable. 

For all 1 1 galaxies in our sample, the choice of host-galaxy 
template has almost no effect on the derived intrinsic SED short- 
ward of A,Brk- On the other hand, in at least five cases the SED at 
longer wavelengths is strongly dependent on the choice of host- 
galaxy template, meaning the exact fom of the intrinsic AGN SED 
at these wavelengths is uncertain. However, in all cases the intrin- 
sic SED falls rapidly (in vFv) beyond 15-60 fi.m, irrespective of 
the choice of host-galaxy component. Despite this fall-off at long 
wavelengths, the total emission at 60 |J.m is still dominated by the 
AGN in at least three, possibly four, cases (namely, MCG-03-34- 
064, ESO-103-035, IC5063 and, possibly, Mrk 3; there is one uns- 
elected host-galaxy component that dominates at 60 \xm. in this last 
case), irrespective of the choice of host-galaxy SED (see fig.|5}. 

In fig.|6]we plot the full range of possible intrinsic AGN SED 
components extracted from our models. Shown in this plot are the 
intrinsic SEDs from all suitable fits identified using the criteria out- 
lined above. A notable feature of this plot is the large range of in- 
trinsic SED shapes extracted from the fits to the 11 observed SEDs 
in our sample. For example, when normalised at 19 \im the flux 
the intrinsic AGN flux extracted from our model fits spans over an 
order of magnitude at 60 \xm. and 100 \im. While the range mid- 
infrared (i.e., A. < 25 (xm) SEDs is well constrained by the AGN- 
dominated IRS spectra, it is not clear how much of the scatter at FIR 
wavelengths (i.e., ^ > 25 iJ.m) is due to differences in the true intrin- 
sic SED and how much is introduced by the different host-galaxy 
templates we use. Therefore, the range shown here will cover a 
larger spread than the true intrinsic AGN SEDs at A > 25 \im. In fig. 
[6] we have discriminated between AGNs with X-ray luminosities 
above and below the median of the sample (i.e., above and below 
log(L2-10kev) =42.9). We identify a clear trend for more intrinsi- 
cally luminous AGNs to have intrinsic infrared SEDs that fall more 
rapidly at longer wavelengths, although stress that there is some 
overlap between the SEDs of AGNs with X-ray luminosities above 
and below the median of our sample (indeed, the most rapidly 
falling SED is that of the lowest luminosity AGN in our sample, 
NGC 21 10). We also find that, on average, the mid-infrared SEDs 
of the more X-ray luminous AGNs in our sample have weaker sil- 
icate absorption and stronger silicate emission, although this could 
be related to the fact that all of the low X-ray luminosity AGNs are 
classed as either Type 1.8, 1.9 or 2 AGNs which, as we have seen, 
tend to have stronger silicate absorption features (see i]4. U . The 
evidence for a link between the shape of the intrinsic infrared SED 



and L2_iokeV is strengthened by differences between the average 
intrinsic SEDs of the moderate luminosity AGNs studied here and 
those of more luminous quasars (e.g., the average quasar SED pre- 
sented in Netzer et al. 2007; see i]5.2b . There is weak evidence of 
a relation between the shape of the SEDs shortward and longward 
of ^Brk> with SEDs with lower values of ai falling more rapidly 
at FIR wavelengths. However, we note that there are at least two 
objects in our sample that deviate strongly from this trend (namely, 
Mrk 3 and NGC 5506; see fig.ig. 

We also plot in fig.|6]the average intrinsic SEDs calculated by 
taking the mean of the intrinsic AGN components extracted from 
the fits. Because of the considerable overlap in shape of the intrin- 
sic SED at MIR wavelengths between Type I and Type 2 AGNs 
(see i]4. H and fig.|2} and the small number of intrinsic SEDs in our 
sample (i.e., 1 1) we do not discriminate between AGN types when 
producing these averages. Although not included in this plot, we 
note that the average of the SEDs obtained by subtracting the host- 
galaxy components from the observed data (i.e., those shown in fig. 
O is consistent with the average SED shown in fig.|6] This average 
SED can be expressed as: 

fX'-^ at 6\im<'k<\9 \im 
X"-^ at 19 ^m < X < 40 ^m (I) 
V^'^F^^ at A,>40nm 

and shows evidence of weak silicate absorption (equivalent to an 
absorbing column of A'h~5x lO^^cm^^, or Xqj ~0.2, using fig. 
10 of Draine 2003 to convert A'h to tg 7) but little or no silicate 
emission. For clarity, we plot the average mean-average SEDs of 
log(L2_iokev) >42.9 and log(L2_iokev) <42.9 AGNs separately 
in the right-hand panel of fig.|7] As expected from the trend identi- 
fied above, the average intrinsic SED of the more luminous AGNs 
in our sample (i.e., log(L2_i()kev) >42.9) falls more rapidly at 
longer wavelengths than that of the lower luminosity AGNs (i.e., 
log(L2_iokev) <42.9). For example, when normalised at 19 \im the 
lower luminosity AGNs emit, on average, 2-3 times more flux at 
60 |j.m and 100 iJ.m than the higher luminosity AGNs in our sam- 
ple. The parameters describing the average SEDs of the high and 
low luminosity AGNs in our sample are largely the same as shown 
in Eqn.[T] although the spectral indices at X, > 19 (xm are somewhat 
different, i.e., a2=0.0 and 0.4, respectively. This relative difference 
between the FIR SEDs of high and low luminosity AGNs in our 
sample could be the result of higher luminosity AGNs being capa- 
ble of heating a larger fraction of their surrounding dust to higher 
temperatures. A consequence of this increased heating would be 
relatively stronger emission at MIR wavelengths; hence the appar- 
ent flux deficit at FIR wavelengths when normalised at 19 |J.m (see 
also ^5.2\ . All three average intrinsic AGN SEDs described here 
are published in columns 2-4 of table [3] (which is available in its 
entirety online at http://sites.google.com/site/decompirl. 

Finally, we explore whether using a more straightforward ap- 
proach to normalise the host galaxy component has any effect on 
the extracted intrinsic AGN SEDs. To e xtract the average in trin- 
sic SED of more luminous, quasar AGNs. lNetzer et al1 ( l2007h nor- 
malised the host galaxy component such that it constitutes the ma- 
jority of the observed flux at 60 |J.m and 100 \im. If we take a similar 
approach and assume that 90% of the flux at 100 |J.m is emitted by 
the host-galaxy, we obtain largely the same overall intrinsic SED 
shapes as shown in fig.|5]and fig.|6l but are left with strong features 
at 7.7 [xm where the PAH features are over-estimated by the host- 
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Figure 5. The intrinsic mid-infrared spectra (light blue Hnes) and far-infrared photometry (light blue points) of the 1 1 AGN-dominated sources in our sam- 
ple after subtracting suitable host-galaxy components. The suitability of the host-galaxy components is assessed from the SED fits described in j|4.2| (see 
also appendix|A)- We highlight the SED produced by subtracting the most suitable host-galaxy component (i.e., the "best-of-the-best"; see i|4.2l dark blue 
line/points). The oiiginal observed data is shown in red. The intrinsic AGN infrared SEDs shown here cover a broad range of shapes, although all fall rapidly at 
wavelengths longwards of 15-60 i^m. Despite this fall-off at long wavelengths, there are at least four observed AGNs SEDs (namely, Mrk 3, MCG-03-34-064, 
ESO-103-035 and IC5063) that are AGN-dominated even at 60 |xm, irrespective of the choice of subtracted host-galaxy template. 



galaxy component, indicating that this approach is less reliable than 
our adopted approach^ 



5 COMPARISON WITH PREVIOUSLY DEFINED AGN 
INFRARED SEDS AND RESULTS FROM DUSTY 
TORUS MODELS 

We have produced a set of intrinsic AGN infrared SEDs by decom- 
posing the observed SEDS of local, AGN-dominated galaxies into 



well-defined host-galaxy and intrinsic AGN components (see ^4.21 . 
The full variety of these intrinsic AGN SEDs is shown in fig.|6] In 
this section we explore how these intrinsic SEDs compare against 
other commonly assumed AGN SEDs, including those produced by 
radiative transfer models of the dust surrounding the AGN. For sim- 
plicity, the comparisons made in this section are against the mean 
average and range of all the intrinsic SEDs in our sample (i.e., we 
do not differentiate between the high and low luminosity AGNs in 
our sample; see previous section). 



We also get broadly similar results if we assume that 80% or 99% of the 
100 |j.m flux is due to the host-galaxy. 
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Figure 7. Left: The mean-average and range of intrinsic AGN infrared SEDs from tliis study plotted against tlie observed infrared SEDs of Mrk 23 1 and 
NGC 1068 (i.e., two galaxies commonly used to represent the typical AGN infrared SED). When all are normalised at 19 ^im, the SEDs of NGC 1068 and 
Mrk 231 lie above the intrinsic AGN SEDs at FIR wavelengths, which we attribute to host-galaxy contamination in the two comparison SEDs (see ^5.1\ . Right. 
The mean average intrin sic SEDs of (a) all 1 1 AGNs in our sample, (b) log(L2-iokev)<4-2.9 AGNs and (c) log( L2-iokev)>42.9 AGNs, plotted against the 
average quasar SEDs of Richards et al. (2006) and the average intrinsic quasar infrared SED of lNetzer et al.lfcOO?!) . Again, all SEDs are norma lised at 19 |i.m. 
Note t hat the Richards et al. (2006) and Netzer et al. (2007) SEDs extend to ~ 95 ]im and ~65 p.m, respectively. At this normalisation, the iRichards et al.l 
i2006l) average quasar SED is well matched to the average intrinsic SE D at X >30 |j.m, although it lies above our average intrinsic SED at shorter wavelengths. 
At ^ <19 |xm the average intrinsic quasar SED of iNetzer et al.l ( |2007|) lies slightly above all the average intrinsic SEDs from this study, but lies below them at 
X >19 |J.m. This plot clearly shows the trend reported in i]4.2l for more luminous AGNs to have bluer intrinsic infrared SEDs. 



5.1 NGC 1068, Mrk 231 



5.2 Infrared Quasar SEDs 



Of all known local AGNs, NGC 1068 and Mrk 231 are com- 
monly used to characterise the infrared SEDs of AGNs. In fig. [7] 
we illustrate how our range of intrinsic AGN SEDs compare with 
these canonical AGNs0 When normalised to 19 vim, the SEDs 
of Mrk 231 and NGC 1068 are, respectively, ~15 and ~8 times 
higher at 100 vim than the average SED calculated from our AGN 
templates. Assuming that the average intrinsic AGN emission con- 
tinues to fall as a modified black body to 1000 yim, the Mrk 231 
and NGC 1068 SEDs contain roughly 1.5-3 times the power at 
infrared wavelengths than the average intrinsic AGN SED (when 
all are normalised at 19 |i.m). This excess power is likely the 
result of host galaxy contamin a tion in these two ca s es, as sug- 
gested by iTelesco et alj ( Il984h : iDownes & SolomoiJ ( Il998h and 
iLe Floc'h et al.l ( l200ll) . Indeed, fits to the infrared SEDs of Mrk 23 1 
and NGC 1068 using the approach outlined in i]4.2| are consistent 
with them being host-galaxy dominated at FIR wavelengths. Fi- 
nally, we note that at wavelengths < 20 [im, where the host-galaxy 
contributions are small, the SEDs of NGC 1068 and Mrk 231 are 
largely consistent with our intrinsic AGN SEDs. 



^ The NGC 1068 and Mrk 231 SEDs shown here represent the infrared 
emission produced by the entire galaxy (i.e., including emission from 
the host galaxy). We note that, whilst mid-infrared fl ux density measure- 
ments of the reso lved AGN core of NGC 1068 (e.g., iGandhi et al.ll20O9l ; 
IPrieto et alj|2010l) are now available, we choose to include emission from 
the host galaxy to allow easier comparison with previous studies which 
have, in general, done the same. 



There have been a number of attempts to constrain the aver- 
age quasar SED, although only a handful have included coverage 
to FIR wavelengt hs (e.g. lElvis etaLlll994l : [Richards et al]|2006l: 
iNetzer et al."2007^. One of the most prominent, recent studies is 
that of Richards et al. ( 2006), which combined photometric data 
covering the radio to X-ray regimes for a sample of optically- 
selected, broad-line quasars. In fig.|7]we compare the average in- 
frared SED of all quasars from that study t o the AGN templates 
derived here. The average quasar SEDs from 'Richards et al.l(l2006l) 
turns over at approximately the same FIR wavelengths as the aver- 
age of our intrinsic AGN SEDs, although we note that host-galaxy 
contamination at FIR wavelengths has not been removed from 
this average quasar SED. This host-galaxy contamination will tend 
to push the po sition of the turn-over to longer wavelengths. The 
[Richards et al.l quasar SED is flatter than the average intrinsic SED 
of more moderate luminosity AGNs, which could be due to the in- 
creased quasar luminosities heating the surrounding dust to higher 
temperatures compared to more moderate l uminosity AGNs. De- 
spite being flatter at short wavelengths, the [Richards et al.l ( |2006|) 
SED contains only '-^lO per cent more flux over 8-1000 yim than 
the average intrinsic SED when both are normalised to 19 [im and 
extrapolated to longer wavelengths assuming a modified blackbody 
SED. 

Recently, INetzer et al.l ( |2007|) removed the host-galaxy com- 
ponent from the average infrared SED of a sample of PG-quasars 
to produce an average, intrinsic quasar infrared SED. This SED 
is also investigated in fig.|7| When both are normalised at 19 [im 
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Figure 6. The full range of possible intrinsic AGN infrared SEDs produced 
by our SED fitting procedure described in ^4.21 We include in this plot 
the intrinsic SEDs extracted from the fits deemed suitable using the cri- 
teria outlined in i|4.2l (i.e.. those highhghted in appendix |A). Each intrin- 
sic SED is normalised at 19 \xm to demonstrate the range of SED shapes. 
There is a clear systematic difference between the intrinsic SEDs of AGNs 
with X-ray luminosities above and below the median of our sample (i.e., 
log(L2-l0kev) =42.9). Also included in this plot is the mean-average intrin- 
sic SEDs of all 1 1 AGNs in our sample, labelled to illustrate the parameter- 
isation outlined in Eqn.[T] 



the average intrinsic quasar SED emits more strongly at X < 19 iJ.m 
than the average intrinsic SED of more moderate luminosity AGNs 
calculated above. However, this situation is reversed at A. > 19 |J.m, 
with the former falling more rapidly at longer wavelengths. This 
strengthens our findings described in ^4.2\ where we report a simi- 
lar trend between the higher (i.e., log(L2_iokev) >42.9) and lower 
(i.e., log(L2_iokev) <42.9) luminosity AGNs within our sample. 
We note that the intrinsic quasar falls even more rapidly than the 
average intrinsic SED of the higher luminosity AGNs in our sam- 
ple. The intrinsic quasar SED has stronger silicate emission fea- 
tures than the average intrinsic SED of more modera t e lum inosity 
AGNs, although this may be due to the lNetzer et al.l ( |2007|) being 
comprised of Type 1 AGNs which tend to have stronger silicate 
emission features (see ^4.1\ . As we suggest in i|4.2| the increased 
flux at short wavelengths and relative deficit of flux at longer wave- 
lengths in the average quasar SED may be related to the higher lu- 
minosities of quasars, with more luminous objects capable of heat- 
ing more dust to higher temperatures. However, we can only claim 
this in an average sense as the average quasar SED is consistent 
with the spread of intrinsic infrared SEDs extracted from the 11 
AGN-dominated sources above. The four intrinsic AGN SEDs that 
bear the closest similarity to the average quasar intrinsic SED are 
3C120, IC4239A, Mrk 509 and NGC 4507; three are type 1 AGNs 
and one is a type 2 AGN, respectively. All four of these SEDs have 
intrinsic 2-10 keV luminosities greater than the median of our sam- 
ple, (i.e., Iog(L2_iokev) >42.9)^ 

For comparison, we fit the lNetzer et al.l ( l2007b intrinsic quasar 
SED using DECOMPIR (see appendix|B]l, obtaining: 
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Figure 8. The range of the intrinsic AGN infrared SEDs (as shown in 
fig. [T] left panel) plotted against the range of SEDs predicted by contin- 
u ous (left-hand plot ) an d clumpy (right-hand plo t) torus models described 
in lFritz et alj j200^ and lSchartmann et al.l fcOOsh . respectively (shaded re- 
gions in the left and right panels, respectively). Each set of SEDs has been 
normalised at 19 |xm. In general, a clumpy torus model provides a better 
representation of the intrinsic AGN SED than the continuous torus models 
which over-predict the range of intrinsic AGN infrared SEDs by at least one 
order of magnitude at ~60 |j.m with this normalisation. However, none of 
the clumpy torus models predict 6-20 |j.m SEDs as steep as some of those 
in our AGN-dominated sample (e.g., Mrk 3) 



,1.2 



at 6 \im <X<20 \im 
at A. > 20 yun 



(2) 



where the symbols are the same as those defined in ^4.2\ These 
parameters confirm the shallower S ED at MIR wavele ngths and 
shorter-wavelength turnover of the iNetzer et al.l (l2007h intrinsic 
quasar SED, compared to the average AGN SEP defined in this 
study. When we extrapolate the lNetzer et al. I( l2007l) mtrinsic quasar 
SED to longer wavelengths as a modified blackbody and nonnalise 
both SEDs at 19 yum, then we find that the total power radiated at 
infrared wavelengths is --^80 per cent that of the average intrinsic 
infrared SED of typical AGNs calculated here. 



5.3 Radiative transfer (dusty torus) models 

Finally, we compare the range of intrinsic, infrared AGN SED tem- 
plates with predictions from radiative transfer models of the dust 
surrounding the active nucleus, which is thought to be the main 
source of infrared emission from AGNs (although there is some 
evidence to suggest that infrared emission from AGNs could also 
be produced in more e xtended, dusty regions; e.g. lSchweitzer et al] 
l2008l ; lMor et al.ll2009l) . We compare our templates with the SEDs 
produced by both (a) continuous and (b) cl umpy distribution s of 
dust (specifically, the m odels described in iFritz et all |2006| and 
ISchartmann et al.ll2008l respectively). The range of input param- 
eters of the form er (i.e., continuous) model are given in table 1 of 
iFritz et al] ( |2006|) . The inp ut parameters of th e clump y torus model 
is presented in table 1 of ISchartmann et si\ ( |2008|) ; we used the 
SEDs produced by varying both orientation angle (see their fig. 4) 
and dust mass (see their fig. 10). The SEDs predicted by these torus 
models are shown in fig. jS] together with the average and range of 
intrinsic SEDs defined in this study. 

The continuous dust-distribution models produce a much 
broader range of SED shapes than the intrinsic AGN infrared SEDs 
defined here. For example, when normalised at 19 [im these mod- 
els predict a range of flux densities spanning almost three orders 
of magnitude at 60 [im, compared to a spread of (at most) two or- 
ders of magnitude for the intrinsic SEDs defined in this study. We 
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Figure 9. The intrinsic 1 2 |j.m AGN luminosities derived from our fits to tlie 
IRAS pliotometiy plotted against the core 12 [xm luminosities obtained us- 
ing high spatial resolution (i.e. , sub-arcsecond) m id-in frared observa t ions o f 
23 nearby AGNs, taken from lHorst eTaD l l2008h and lGandhi etal.l j2009h . 
In the latter, the host galaxy contribution is spatially resolved, providing 
an uncontaminated measure of the intrinsic AGN infrared luminosity at 
12 |j.m. After fitting the large-aperture (i.e., arcminute) IRAS photometry 
and excluding any host galaxy contamination we reproduce, on average, 
the intrinsic luminosities derived from high spatial resolution observation 
to within a factor of ~2, compared to a factor of ~4 when the 12 |a.m IRAS 
photometry is not corrected for host-galaxy contamination. Points associ- 
ated with Compton-thick AGNs have been marked with a "C" and show a 
similar level of scatter as the Compton thin sources. 



note that the model SEDs that match the range of intrinsic AGN 
SEDs do not correspond to a particular region of input parame- 
ter space (i.e., torus opening angle, optical depths, inner and outer 
torus radii, radial dust density distribution, angular dust density dis- 
tribution, SED of incident radiation). Therefore, our intrinsic AGN 
SEDs are unable to directly constrain acceptable input parameters 
of these continuous dust-distribution models. 

The SEDs produced by the clumpy torus models provide a 
closer match to the observed intrinsic AGN infrared SEDs, al- 
though they too over-predict the range of FIR fluxes spanned by 
our intrinsic SEDs and tend to peak (in vFy) at longer wavelengths. 
The closer match to the clumpy torus models is consistent with 
other evidence in s u pport for such models (e.g., iHorst et al.ir2006l : 



ottier evidence m s up| 
lGandhietal]|2009l :[F 



Ramos Almeida etal]|2009t) However, none 



of the clumpy torus models described in lSchartmann et al] dlOOSl) 
predict 6-20 (i.m SEDs as steep as some of those in our AGN- 
dominated sample (e.g., Mrk 3). As a consequence, at this normal- 
isation the clumpy torus models over-predict the emission of many 
of our sources at X < 19 iJ.m. 



6 TESTING THE AGN AND HOST-GALAXY 
TEMPLATES 

A principal motivation behind defining the intrinsic infrared SEDs 
of moderate luminosity AGNs is to use these as templates to fit the 
infrared SEDs of composite galaxies and calculate the AGN and 
host galaxy contributions to their total infrared output. This is espe- 
cially useful in the majority of cases when only broad-band infrared 
photometry measurements are available to constrain the infrared 
SED. In this section we explore whether fits to the broad-band pho- 
tometry incorporating the host-galaxy and AGN templates defined 
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Figure 10. EW PAH^^n 25 plotted against the AGN contribution at 19 |a.m 
derived from our fit s to IRAS photometry of 78 sources in the 12 |x m sample 
of Seyfert galaxies jRush et al.l 199llTommasin et alj2008ll2010ll . Sources 
where the IRAS Yl ycm. photometry measurement is at least 2 times higher 
than the flux at the same wavelength in the high resolution IRS spectrum 
in which the EW PAHj^n 25 is measured have been indicated. Such cases 
are likely to suffer from significant aperture effects which will typically 
lead to more of the host galaxy flux being included in the IRAS photome- 
try measurements. Also shown in this plot is the EW PAH^ 1 1 9^ - percent 
AGN contribution relation taken from iTorrmiasin et al.l fcOlOb (black line) 
with a ±25 per cent margin shown in grey. We note a general agreement 
between these two measures of the host-galaxy contribution, with ~ 76 per 
cent (i.e., 59/78) of AGN contributions measured from the photometry fits 
lying within 25 per cent of that expected using EW PAHj^n 25. 



here can reproduce the intrinsic AGN luminosities and AGN contri- 
butions derived using other, independent approaches (i.e., high res- 
olution imaging of nearby AGNs and emission line diagnostics). 
To perfoiTn these fits, we have developed an IDL procedure, DE- 
COMPIR, which we describe in appendix|B]and is publicly available 
online at |http://sites.google.com/site/decompir| 



6.1 Comparison with high spatial resolution observations of 
local Seyferts. 

A number of recent studies (e. g. iKrabbe et al.ll200ll ; iHorst et all 
I2OO6I I2OO8I; iGandhi et alj|2009l) have shown that the intrinsic in- 
frared luminosities of nearby AGNs can be reliably measured from 
high spatial resolution observations (i.e., ~ 0.5", typically probing 
physical scales of < 100 pc). At such high resolutions, host galaxy 
contamination is minimised, leaving only the emission intrinsic to 
the active nucleus. However, such analyses are only viable for a 
small number of systems that are nearby enough that the central 
regions can be sufficiently well resolved. On the other hand, if we 
can use SED decomposition to exclude the host galaxy contribu- 
tion to the broad-band SEDs, then we could potentially measure 
the intrinsic 12 (j.m luminosities for all AGNs with well sampled 
infrared SEDs. Here, we demonstrate that the SED decomposition 



procedure outlined in this study can, indeed, be used to measure the 
intrinsic AGN 12 \im flux. 

We use DECOMPIR (see appendix|Bl( to fit the four-band IRAS 
p hotometries o f 23 A G Ns with core 1 2 ^.m luminosities reported 
in lHorst et iD ( l2008[) or lGandhietai] ( |2009|) and X-ray luminosi- 
ties within the range spanned by the AGNs used to define our in- 
trinsic AGN SEDs (i.e., L2_iokeV= lO'*^ - 10'*'* ergs s"*). As we 
are limited to only four independent flux density measurements we 
only allow the normalisations of the AGN and host galaxy compo- 
nents to vary when fitting the observed SEDs. All other parameters 
are fixed to their average values derived from the full sample of 
11 AGN dominated galaxies described above (also see appendix 
iBt . but we note that our main results do not change if we use the 
average SED of either the high (i.e., log(L2_iokev)>42.9) or low 
(i.e., log(L2-i0kev)<42.9) X-ray luminosity AGNs in our sample. 
Each SED was fit five times, once for each of our host-galaxy tem- 
plates (i.e., as recommended in appendix |B). To calculate the in- 
trinsic 12 \im AGN flux, we integrate the intrinsic AGN compo- 
nent derived from the best fitting solution of these five fitted SEDs 
(i.e., that with the lowest associated value) over a narrow (i.e., 
1 fi.m) top-hat passband centred at 12 [im, which is a similar re- 
sp onse function as th ose filters used for the ob servations described 
in lHorst et all ilOO^) and lGandhi eTai] ( l2009h . 

In fig. [9] we plot the intrinsic 12 [im flux derived from our best 
fitting SEDs against those obtained usin g high spat ial resolution 
obser vations as reported in iHorst et ^ |2008) and iGandhi et al.] 
( I2OO9,) . We find that the intrinsic luminosities derived using our 
SED fits are well matched to those obtained using high resolution 
observations of the AGN core. On average, the intrinsic 12 (im 
AGN luminosity derived from our fits lie within a factor of two 
of the core luminosity at this wavelength, compared to within a 
factor of four when the 12 [xm IRAS photometry is not corrected 
for host-galaxy contamination. A linear regression to these points 
gives: 



rFil 
^12 



erg s 



= (0.13±0.07)- 



(0.96±0.09)log(^^.^^_ 



(3) 



which has a slope consistent with a 1:1 relationship. For compari- 
son, we also plot the 12 |j.m luminosities calculated directly from 
the 12 p.m IRAS photometry, noting that these often overestimate 
the intrinsic 12 [im luminosities, sometimes by more than a factor 
of 5. Finally, we note that the five Compton-thick sources in this 
sample (i.e., N GC 1068, IC 363 9, Swift J0601. 9-8636, NGC 5728 
and NGC3281: lDella Ceca et a l. 2008) display a similar degree of 
scatter as the full sample. This demonstrates the power of using 
infrared SED fitting to identify even heavily obscured AGNs. 



used for this section of our study, together with their measured 
emission line fluxes). Again, as we are dealing with only four pho- 
tometry measurements, we only allow the normalisations of the 
AGN and host-galaxy components to vary when fitting the ob- 
served data. All other parameters are fixed to their average values 
derived from the full sample of 1 1 AGN-dominated galaxies (see 
i]4.2l again, we note that our results do not change if we use the av- 
erage SEDs derived from the high (i.e., log(L2_i()kev)>42.9) or low 
(i.e., log(L2-i0kev)<42.9) X-ray luminosity AGNs in our sample). 
The IRAS photometry of each of the 78 objects in our sample are 
fitted five times, each time using a different choice of host-galaxy 
template (defined in i|3j. Here, we report the AGN contribution to 
the 19 p.m infrared emission measured from the best fitting of these 
five independent fits (i.e., that with the lowest associated value). 

In fig. [To] we plot the AGN contribution at 19 ^im derived 
from our fits to the IRAS photometry agai nst E^V PAH) ^ 1 mea- 
sured from high resolution IRS spectra; seel Tommasin et aI1 ( l2008l . 
I2OIOI) . We find that the AGN contribution measured from our fits 
to the IRAS photometry are in good general agreement with those 
based on EW PAH;j^ii.25' especially at high levels of AGN contribu- 
tion (i.e., > 50 percent). Of the entire sample of 78 AGNs, ~76 per 
cent (i.e., 59/78) have AGN contributions measured from our broad 
band SED fits that lie within ±25 per cent of those derived from 
EW PAH;,^!! 25- At lower AGN contributions where the intrinsic 
AGN component is more difficult to constrain, our fits typically es- 
timate the AGN contribution to within ~ 50 per cent of that derived 
from EW PAH^j^n 25- We note, however, that at such low AGN con- 
tributions even emission line diagnostics have considerable diffi- 
culty in estimating the intrinsic AGN contribution, as demonstrated 
by the significant amounts of scatter in the [Ne v]/[Ne S3. 2.2 II 
- EW PAHj^i] 25 coiTelation shown in fig. 4 of Tommasin et al 
(120 lOT) (see also right-hand panel of fig. 8., iGoulding & Alexandei 
l2009t) 



7 APPLICATION OF OUR TEMPLATES AND 

INFRARED PHOTOMETRY FITTING PROCEDURE 

In the previous section, we quantified the accuracy of our SED 
fitting routine by comparing the results from these fits to results 
obtained by other well-established approaches In this section, we 
demonstrate two additional applications of our AGN intrinsic AGN 
SEDs and fitting approach: (a) defining the correction factors to 
convert 12 [xm and 2-10 keV luminosities to total AGN infrared 
luminosities (i.e., i^i^''^) and (b) measuring the intrinsic AGN lu- 
minosities of large samples of composite galaxies for which only 
infrared photometry measurements are available. 



6.2 Comparison with results derived from emission line 
diagnostics 

A number of infrared emission line diagnostics have recently been 
identified as providing reliable proxy measures of the AGN and 
host-galaxy contributions to the infrared and bolometric output 
of co mposi t e galaxies (e.g. [G enzel et al] 1 19981 : iTommasin et al.l 
l2008ll20I(]| : [Goulding & Alexa nder 2009). A key test for our fits 
is whether they can broadly reproduce the intrinsic AGN luminosi- 
ties derived from these diagnostics. 

We use DECOMPIR to fit the IRAS photometries of 78 galax- 
ies from the IRAS 12 [im sample for which accurate emission line 
flux es, measured from high reso lution IRS spectra, are available 
(see lTommasin et alj200al2010l for a complete list of all the AGNs 



7.1 vZ^(12 |xm):L^'^ and L2_iokev:^iR''^ correction factors 

In cases where the infrared SED is either poorly constrained or 
dominated by emission from the host galaxy, measuring the total 
intrinsic AGN contribution can prove extremely difficult. Further- 
more, to accurately estimate key parameters such as star-formation 
rates from infrared emission requires that any AGN contribution be 
accounted for. However, often, the only accurate way to measure 
the levels of AGN activity in composite galaxies is via X-ray or 
high resolution infrared observations. Using the intrinsic AGN in- 
frared SEDs defined here coupled with the correlation between the 
in trinsic 12 \im lumin osity and the 2-10 keV luminosity defined 
in lGandhi et al. II2OOS), we can now define a coiTection factor to 



Table 3. Average AGN and host-galaxy SEDs and the wavelength-dependent correction factors used 
to calculated LJ^p^^and Lli^°. 
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Notes: (1) Wavelength, (2)-(4) Flux densities of the average intrinsic AGN infrared SEDs defined 
in [|4] normalised at 6 [xm. Hi. Lum and Lo. Lum refer to the average SEDs of log(L2-lOkcv)> 42.9 
and log(L2-iokev)< 42.9 AGNs in our sample, respectively. (5)-(9) Flux densities of the five host 
galaxy templates defined in ij3] normalised at 6 |xm, (10) Conversion factors to convert vL^™ to total 
intrinsic AGN infrared luminosities (i.e., ij^^) at all wavelengths covered by our templates, based 
on the averag e AGN SED def ined in [|4] (11) The same as column 10, but for the average quasar 
template of Netzer et al. (2007). This table is available in its entirety at higher wavelength resolutions 
at http ://sites/google . com/ site/decompir 



convert 2-10 keV luminosities to total intrinsic AGN infrared lu- 
minosities. This correction factor enables any AGN contribution to 
the infrared output of galaxies for which, for example, only X-rays 
or monochromatic infrared observations provide a measure of the 
intrinsic AGN luminosity. 

In fig.[TT]we show the range of L^^/vZ^( 12 \xm) ratios for 
all 11 AGNs for which we measure the intrinsic AGN infrared 
SED, plotted against vLv(12 nm). Here, we have assumed that the 
intrinsic AGN infrared SED continues as a modified black body 
to 1000 [im. These corrections therefore only apply to the ther- 
mal component of the intrinsic AGN SED as we do not consider 
any emission from non-thermal (e.g., synchrotron) processes. We 
choose vLv(12 [im) as a reference point as it has been shown to 
be st rongly correlated w ith the intrinsic luminosity of the AGN 
(e.g., iGandhietal .1120091) . However, as we have defined the intrin- 
sic AGN SED at 6-1000 [im, an equivalent figure/analyses could 
be made/performed using any wavelength in this range as a refer- 
ence point (see tableO. The L^^/vLv( 12 p.m) ratios displayed by 
these 1 1 AGNs span the range ~ 1.3 to ~3.5, with an average ratio 
of ~ 2.2 (derived from the average SED defined in gj. We do not 
see any correlation between L^^/vLy (12 )i.m) and vLy ( 1 2 l-tm) for 
our sampl e. The averag e Lf^^/vLy{l2 ]mi) ratio for quasar SEDs 
defined in iNetzer et al.l ( [2607i) is slightly lower (i.e., ~ 1.6) than 
that of the lower luminosity AGNs considered here, although this 
difference is not significant based on the range of ratios spanned by 
our sample of intrinsic AGN SEDs . How ever, we note that the at 
other wavelengths the iNetzer et al.l ( l2007ft quasar bolometric cor- 
rection differs considerably from the average bolometric correction 
of more modest luminosity AGNs defined here (see table |3}. We 
also show in fig.[TT]the Lir/vLv(12 (xm) ratios for NGC 1068 and 



Mrk 231, the latter of which is significantly higher than the range 
of intrinsic ratios due to the additional host-galaxy components in 
these SEDs (see ^5H . 

Using the vL v(12 jxm):L2-i0keV correlation defined in 
iGandhi et Zl ( |2009|) , we convert our L^^/vLv( 12 ^.m) ratios into 
L^'^/L2-i0keV correction factors. Assuming L^'^/vLv(12 Vim) = 
2. 2 ± 1.3 (to cover the full range of ratios displayed by our sample) 
we obtain: 

/ T AGN \ 

l°g(T5^)= (0-53 ±0.26)+ 

(l.ll±0.07)log(jia^) 

For a L2-l0keV= 10^"^ srgs s^' AGN, typical of the range of 
^2-lOkeV considered here, equation |4] gives L2_iokeV^^iR''^~0-3- 
This compares t o a median ratio of L2-i(n^e\//L^^>0.04 for the 
quasar sample of lElvis et al . ( 1994), although in this case there will 
be some contribution to from the host gal axy that has not 

been accounted for, hence the lower limit (see also lAlexander et aH 
l2005l : we have assumed Lir~ 1.3Lfir)- We attribute the factor of 
<7 difference between these ratios to both the unaccounted-for host 
galaxy contribution and the different fractions of the total bolomet- 
ric luminosities of AGNs (i.e., Lboi) that are e mitted at X-ray ener- 
gies b y quasars and more modest AGNs (see IVasudevan & Fabi^ 
l2007h . For example, i2-l0keV= 10'*^ ergs s ' AGNs emit approx- 
imately 5 per cent of their bolometric power in the 2-10 keV band, 
compared to just ^ 1 per cent for a L2_iokeV= 10^^ ^^S^ ^ ' 
quasar (see fig. 3 of Vasudevan & Fabian 2007 |; i.e., changes in 
apx with lununosit y;[Avni & Tananba um.l98Z : iWilkes et aljl994 
IVignali et alj l2003l : [steffen et al. 200^- 

Using these conversion 
factors, we find that Lir/LboI are similar for both moderate and high 
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Figure 11. The range of L^^/vLv( 12 [im) ratios plotted against 
vLv(12 |j.m) for the 11 AGNs that we use to define the range of intrinsic 
AGN infrared SEDs (filled circles). We also show the L^™/vLv(12 [im) 
ratio for the average intrinsic AGN infrared SEP define d in [|4] the av- 
erage intrinsic quasar infrared SED from Netzer et al.l (2007) and the 
two arch etypal AGNs, NG C 1068 an d Mrk 231 (the average quasar 
SEDs of lElvisetal]|l994l and [Richards et al. 2006 have approximately 
the same Lf^^/vLv{12 |j.m) ratios as the Netzer et al.. .2Q07 quasar SED; 
i.e., L;^™/vLv(12 ^m) «1.60, 1.59 and 1.62, respectively). On av- 
erage, the moderate luminosity AGNs considered in this study have 
higher intrinsic Lj^^/vLv(12 |j.m) ratios than more luminous AGNs (i.e.. 
quasars), although there is no clear correlation between vLv(12 |xm) and 
L^™/vLv(12 [im) for the 11 AGNs in our sample. We wish to point out 
that, as we have defined the average intrinsic AGN SED at all infrared 
wavelengths (i.e., 6-1000 ^m), an equivalent plot could be produced for 
any wavelength in this range (see table|3j 



luminosity (i.e., quasar) AGNs (i.e., 0.17 and <0.25, respectfully). 


Finally, as we have defined the average intrinsic AGN in- 
frared SED at 6 \im to 100 ^m, we can use it to convert an in- 
trinsic AGN infrared flux obtained at any wavelength covered by 
our SED templates into i^'^- We publish these conversion factors 
in column 8 of table [5] (which is available in its entirety online at 
[http:// sites . google.com/ site/decompir) . 



7.2 Application of the SED fitting procedure to the 
Swift-BAT AGN sample 

One of the principal applications of our SED fitting procedure is 
to measure the AGN contribution to the infrared output of compos- 
ite galaxies. Because our fitting procedure can work with sparsely 



* If we assume that Lir/Lbqi provides a measure of the covering fraction 
of the infrared emitting dust, then these numbers are at odds with other 
measures based on the relative numbers of Type 1 and Typ e 2 AGNs in the 
local Universe (i.e., 0.5-0.8; e.g., [ Huchra & Burg 1992; Rush et al.ll 19931 ; 
iMaiolino & Riekdl 199 A I Maia et al. 2003). However, clumpy torus mod- 
els predict that many additional parameters (such as the temperature of the 
clouds, the inner and outer radii of the torus and the total number and distri- 
bution of dust clouds; iNenkova et ai]|2002h on top of the covering fraction 
play a role in de tennining Lir/LboI- Indeed, the clumpy torus models of 
ISchartmann et al . (2008) predict a range of Lir/LboI spanning ~ 0.15 — 0.3 
despite their models all having a constant covering fraction of 0.5. 



sampled photometry data (e.g. four-band IRAS photometry) this 
opens up the possibility of measuring the intrinsic infrared emis- 
sion of AGNs in large samples of galaxies. Furthermore, as the in- 
trinsic AGN luminosity at 12 \Lm is strongly correlated to i-a-lokeV! 
we can use these fits to estimate the intrinsic X-ray luminosities of 
those AGNs that are too heavily obscured to be detected by the cur- 
rent suite of X-ray observatories (e.g., Chandra, XMM-Newton and 
Siizakii). This is analogous to using high-resolution 12 \im obser- 
vations to measure the intrinsic AGN luminosities. To demonstrate 
this, we use the procedure described in appendix |B]to fit the 12 yim, 
25 \im, 60 [xm and 100 ^.m photometry of all of the X-ray AGNs 
in the Swift-BAT sample that were detected in all four IRAS wave- 
bands. 

We use DECOMPIR (see appendix |B]( to fit the IRAS pho- 
tometry of the 44 5vviy!'-BAT AGNs detected in all four infrared 
bands. As we are limited to only four independent flux density 
measurements we only allow the normalisations of the AGN and 
host galaxy components to vary. We fix all of the other parameters 
to their average values derived from the full sample of 1 1 AGN- 
dominated galaxies (see i]4.2| and appendix|B](, but note that the re- 
sults do not change if we instead use the average SEDs of the high 
(i.e., log(L2_iokev)>42.9) or low (i.e., log(L2_iokev)<42.9) X-ray 
luminosity AGNs in that sample. Each SED was fit five times, once 
for each of our host-galaxy templates (i.e., as recommended in ap- 
pendix |BJ. We use fit with the smallest associated value to de- 
termine the AGN contribution to the infrared SED. 

In fig. [12] we plot the intrinsic AGN 12 \Lm luminosities of 
the 42 5w(/f-BAT AGNs against their 14-195 keV X-ray luminosi- 
ties (obtained from lTueller et al.ll2008l) . Here, we use 14-195 keV 
luminosities, rather than 2-10 keV luminosities, as these energies 
are less susceptible to absorption and, as such, provide a more di- 
rect measure of the intrinsic AGN luminosity. Each intrinsic AGN 
12 |j.m luminosity is calculated by integrating the AGN component 
obtained from the best fitting SED over a narrow (i.e., 1 yum) top- 
hat passband centred on 12 yum. We find that the intrinsic 12 yum 
luminosity is strongly correlated with the 14-195 keV luminos- 
ity. For c omparison, we also include in fig. [12] the 42 AGNs pre- 
sented in iGandhi et alj ( |2009|) that have published intrinsic 12 yum 
luminosities measured from high resolution MIR observations. For 
these, we have converted their intrinsic 2-10 keV luminosities to 
14-195 keV luminosities using a constant conversion factor of 2.5 
(i.e., Li4_i95 |jeV=2-5L2_iokev)- We note that there is an inherent 
uncertainty in this conversion factor meaning the true scatter of 
these Li4_i95 j^ev values will differ from that shown in fig. [12] al- 
though we expect any such differences to be comparatively smallQ 
The results derived fr om our SED fits s how a similar level of scatter 
as those presented in Gandhi et al. I( l2009i) . We note, however, that at 
low X-ray luminosities (i.e., < 10^^ ergs s^'), our procedure tends 
to overestimate the intrinsic AGN luminosities compared to those 
obtained from high spatial resolution observations. This is because 
the host galaxy tends to dominate over the AGN at lower luminosi- 
ties, making it increasingly difficult to reliably constrain the AGN 
component. Therefore caution should be taken when applying this 
kind of analyses to studies of lower-luminosity AGNs. A linear re- 



Assuming an X-ray spectral index of r=1.8, we get a factor of 2.7, for 
r=1.3 the factor is 10 and for r=2.3 it is 0.8, implying a factor of ~ 3 un- 
certainty for a typical range of F, although we note that other X-ray compo- 
nents (e.g. cold reflection, high energy turn-over) will also influence these 
factors. 
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Figure 12. vLv(12 |j.m) derived from our DCMB fits to the IRAS pho- 
tometry of 44 AGNs from the Swj/i'-BAT sample plotted against the 14- 
195 keV X-ray luminosities of these AGNs. We see a clear correlation be- 
tween these two measures of the intrinsic AGN luminosities. To give an 
impression of the intrinsic scatter in this relation, we also include 42 AGNs 
which have vLv(12 |xm) measured from high resolution observations. For 
this latter sample, we convert L2-l0kcV to ^14-195 keV using the conversion 
-^-14-195 kcV=2.5L2-10keV- 



gression to all Swift-BAT AGNs considered here gives only modest 
errors, i.e.: 



log 



vL^™(12 |4m) 
10«erg 



(0.37 ±0.08) + 

(0.74±0.13)log 



10«e 



(5) 



with a standard error (i.e., 1 o scatter) of ^ 0.5 dex (i.e., a factor of 
3.2), compared to a standard erro r of 0.36 dex for the 42 com- 
parison AGNs from lGandhietal.l ( l2009l) based on high-resolution 
12 |j.m AGN luminosities. 

The ability of our fitting procedure to derive the intrinsic lu- 
minosities of AGNs from only sparsely sampled photometric in- 
frared data opens up the possibility of using such analyses to ex- 
plore large populations of AGNs detected in wavebands probed by 
the IRAS, ISO, AKARI, Spitzer and Herschel telescopes. Such stud- 
ies will complement work at other wavelengths in exploring the 
properties of both obscured and unobscured AGNs at all redshifts 
probed by these facilities. Applying our SED fitting procedure to 
Spitzer and Herschel data will be the focus of a future study. 



8 SUMMARY 

Using a carefully selected sample of X-ray AGNs that show lit- 
tle or no evidence of host-galaxy contamination in their Spitzer- 
IRS mid-infrared spectra, we have defined the range of intrinsic in- 
frared (i.e., 6-100 \im) SEDs of typical, moderate luminosity (i.e.. 



n42_ 



n44, 



(i.e., z < 0.1). These intrinsic AGN SEDs can be used in conjunc- 
tion with host-galaxy templates to reproduce at least 85 per cent of 
the 6-35 [im IRS spectra of composite galaxies to within ±10 per 
cent as well as their 60 |i.m and 100 [im photometries (see i]4.2l and 
fig. lAlb . We outline a procedure that uses the average intrinsic AGN 
infrared SED, in conjunction with a set of host-galaxy templates, 
to fit the broad-band infrared photometry of composite galaxies. To 



test how accurately this procedure can measure (a) intrinsic AGN 
luminosities and (b) the AGN contribution to the total infrared out- 
put of composite galaxies, we use it to fit the IRAS photometry of 
two independent sam ple of galax i es; th e sa mple of c ore-resolved 
Seyferts described in iHorstetall ( l2008h and lOandhi e t al. (200^ 
and a subsample of the 12 yi m sample of Seyferts (e.g.. , Rush et al.l 
ll993l : lTommasin et al.ll2010h . The results from these tests compare 
favourably against other measures based on high-spatial resolution 
mid-infrared observations and infrared emission line diagnostics. 

We summarise the main results of this study in the following 
points: 



• The intrinsic SED of moderate luminosity (i.e., L2_iokeV = 
1042 _ io44 ergs s^i) AGNs at 6-(15-60) urn (depending on 
turnover point, 15 < A.BB/l-irn<60) can be described as either a sin- 
gle or a broken power-law. When we see a break, which typically 
lies between 15 and 20 |i.m (mean: 19 (i.m), the power-law spec- 
tral indices shortward of the break span the range: 0.7 < ai < 
2.7 (mean: 1.7), where Fy = X". Beyond this break, the intrin- 
sic AGN SED always flattens, with spectral indices spanning the 
range < a2 < 1-5 (mean: 0.7). These ranges of spectral indices 
and and break positions are consis tent with previous studies of the 
mid-i nfrared SEDs of AGNs (e.g.. [Buchanan et al.ll2006l;IWu et all 
l2009l) . See 

• All of the 1 1 AGN dominated infrared SEDs that we can re- 
liably extrapolate to 100 (xm are consistent with being composed 
of a host-galaxy component plus an intrinsic AGN component. All 
of the intrinsic AGN SEDs peak (in vFy) between 15 and 60 \im 
before falling rapidly at longer wavelengths. Despite this fall-off at 
long wavelengths, there are at least 3 AGNs in our sample whose 
observed infrared SEDs is dominated by the AGN even at 60 \un. 
Seeig^l 

• We find evidence that the shape of the intrinsic AGN infrared 
SEDs is related to AGN luminosity, with more luminous AGNs 
having bluer intrinsic SEDs (i.e., more luminous AGNs emit, rela- 
tively, more strongly at MIR wavelengths and less strongly at FIR 
wavelengths) . This trend is seen both within our sample (see ^4.2\ 
fig.[6]and fig.|7]( and when we compare our results with the intrinsic 
SEDs of more luminous quasars (see ^5.2\ . 

• The range of intrinsic AGN infrared SEDs reported here are 
more consistent with the results of clumpy, rather than continuous, 
dusty torus mode ls. These latter models, such as those described in 
iFritz et al. liooi), cover a much broader range of SED shapes than 
the intrinsic AGN infrared SEDs defined here. However, we also 
note that even clumpy torus models have a tendency to over-predict 
the intrinsic infrared emission of AGNs at A. > 30 vim. Furthermore, 
clumpy torus models do not produce 6-20 \im as steep as that of a 
number of AGNs in our sample. See H5.3\ 

• Using the intrinsic AGN infrared SEDs defined in S|4]we de- 
fine a set of correction factors to convert intrinsic AGN infrared 
or X-ray luminosities to total infrared luminosities (i.e., integrated 
8-1000 |j.m luminosities). These conversion factors can be used to 
estimate the AGN contribution to the total infrared output of galax- 
ies when only a single (intrinsic) infrared or X-ray constraint is 
available. See i]7.1l 

• Fits to the infrared SEDs of composite galaxies can be used 
to measure the intrinsic 12 \im luminosities of AGNs to within a 
factor of ~2 when only low-resolution, broad band photometry is 
available ( i]6.1| l. The same fits can be used to determine the intrinsic 
2-10 keV X-ray luminosities to within a factor of ~3 ( i]7.2| ) and 
the fractional AGN contribution to the overall infrared output of a 



composite galaxy to within 25-50 per cent (depending on whether 
the galaxy is AGN dominated; %.2i . 
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APPENDIX A: PLOTS OF FITS TO THE IRS SPECTRA 
AND IRAS PHOTOMETRY 

In fig. lAll we present the observed SEDs of our sample of 1 1 AGNs 
that we use to defined the intrinsic infrared AGN SED. We have 
included in each of these plots the fits to these SEDs produced us- 
ing our five host-galaxy templates (see ^ and an intrinsic AGN 
component (described in i]4.2t . Embedded in each plot is a panel 
showing a zoomed-in image of the fit to the 6-35 |i.m IRS spectrum 
and the residuals produced by subtracting the fit from the observed 
data. 



APPENDIX B: DESCRIPTION OF THE FITTING 
ROUTINE: DECOMPIR 

The results presented in fj6] and ijT] were obtained using a new 
infrared SED fitting routine, DECOMPIR, that we have developed 
which uses the host-galaxy and intrinsic AGN SEDs described in 
S|4]to measure the contribution to the infrared output of composite 
galaxies from these two components. In this appendix, we provide 
a description of the main features of DECOMPIR routine and some 
brief instructions for its use.0 

DECOMPIR is written in the IDL programming language and 
uses minimisation to identify the best fitting model SED to the 
observed infrared data (either spectra and/or photometry) by vary- 
ing the values of a set of free parameters (e.g., component normali- 
sations, AGN component shape, level of dust extinction). DECOM- 
PIR makes heavy use of the IDL package MPFIT to perform this 
fit0 In its default mode, DECOMPIR will attempt to fit the supplied 
infrared photometry or spectrum with the average intrinsic AGN in- 
frared SED defined in i}4]plus the "SBl" host galaxy template (see 
Sj3](. By default, only the normalisations of these components are 
allowed to vary. However, DECOMPIR can be instructed to allow 
any of the parameters used to describe the AGN component to vary 
within the range of values outlined in [j4](i.e., ai, a2, ^biIo ^bb)- 
DECOMPIR can also be instructed to allow the levels of extinction 
applied to both the host-galaxy and AGN components to vary. We 
note that DECOMPIR will only use one host-galaxy template at a 
time (i.e., an SED cannot be fit with multiple host-galaxy compo- 
nents). We therefore recommend that each observed SED is fit five 



° The IDL source code for DECOMPIR has been made available at 
http://sites.google.com/site/decompir 

MPFIT is also publicly available from http://purl.com/net/mpfit [ and is de- 
scribed in lMarkwardll | |2009|) . 
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Figure Al. On this and the following 10 pages we present the 11 AGN-dominated SEDs that we use to extend the intrinsic AGN infrared SEDs to 100 \mi. 
The five panels are used to show the results of fitting the observed SEDs with models containing the five separate host-galaxy components described in ij3] 
In each panel the IRS spectrum (±5o) is shown as grey shading. In each case, the model SED fit is represented as a solid black line. The AGN component is 
shown as a dashed line and the SB component is shown as dotted fine. Where present, the silicate emission at ~10 (j-m and ~18 |^m is shown as a dash-dot 
line on top of the AGN component (see i|4.2l for further details on the fitting procedure). A zoomed-in image of the fit to the 6-35 |J.m portion of the IRS data 
is shown inset, together with the residuals of the fit (the positions of the PAH features at 7.7 \mi, 11.25 i^m and 12.8 |j.m are highlighted with vertical dotted 
lines). We use IRAS 60 \mi and 100 |j.m photometry to constrain the emission at far-infrared wavelengths. We have placed an asterisk next to the starburst 
name of each fit that is considered good based on the criteria outlined in i|4.2| 
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Figure Al - continued 



times, once for each of the host-galaxy templates defined in ^ The 
best fitting SED can then be selected from the five separate fits. 

Because DECOMPIR relies on minimisation to derive the 
best fitting parameters, the numbers of parameters that are allowed 
to vary depends on the number of independent flux density mea- 
surements being fit. Specifically, the number of variable parameters 



must be less than the number of independent flux measurements 
to prevent the fits from becoming degenerate. Some judgement 
based on the nature of the research and the number of independent 
flux measurements is needed in order to decide which parameters 
should be allowed to vary. As a guide, we provide a list of variable 
parameters in order of importance, with the first having the largest 




effect on the fitted SED: tti, dust extinction (botfi components), 

a2, silicate emission features, Xbb, ^Brk (we assume that the rela- 
tive normalisations of the components are always allowed to vary). 
For example, if only five or six independent flux measurements are 
available (e.g. Spitzer 8 (j.m, 24 [xmplus Herschel 100 |im, 160 |J,m, 
250 |im photometry) we suggest that only ai , together with the 



normalisations of the AGN and host-galaxy templates, is allowed 

to vary. We note throughout the main text that the choice of av- 
erage SED template makes little difference to the results from de- 
composing the photometric SEDs of samples of composite galax- 
ies. The reason why this is the case is that, at FIR wavelengths 
(i.e.. A, > 30 \im) the SED is typically (although not exclusively. 
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Figure Al - continued 



see i]4.2t dominated by the host-galaxy. Therefore, a difference of 
a factor of 2-3 between the intrinsic AGN SED at 100 (im has little 
effect on our results when the observed SED is > 80 per cent dom- 
inated by the host galaxy at this wavelength. However, extra care 
must be taken in selecting the appropriate intrinsic AGN SED when 
it is possible that the FIR SED is dominated by AGN emission (i.e., 



when the observed SED falls strongly longwards of ~30 |i.m, e.g., 
MCG-03-34-064, IC5063, ESO103-035, Mrk 3, from this study; 
see fig-O 




10 15 20 25 30 35 
_i 1 1 1 1 1 

10 100 
Wovelength {fim) 



Figure Al - continued 



10" 



10' 



10'' 



10" 



10" 



10" 



10'- 



10'' 



10" 



10" 
10" 



-■ — i i i I 1 

MCG-03-34 -064 



SB1 



5 

» 4 
2 3 

=> 1 


0.2 
0.0 
-0.2 



10 15 20 25 30 35 

H ^ ^ I 

SB3« 





10 15 20 25 30 35 

H — I I I I 1 1 1 1 1 — I I I I 

SB5. 



10'- 



10'' 



10" 



10" 



' J ' ^ s H _ , _ 



5 

» 4 
2 3 

=> 1 

0.2 - 
0.0 
-0.2 




'-I - A 



'MAM 



10 15 20 25 30 35 
_l 1 1 1 1 1 



10" 



10'- 



10'' 



10' 



10" 
10" 



10'- 



10'' 



10" 



10" 



SB2« 




V. 



5 

■o 4 

2 3 

1 


0.2 
0.0 
-0.2 





1 

















10 15 20 25 30 35 



SB4. 





10 15 20 25 30 35 



10 



Wovelength (jim) 



100 



10 



Wovelength (p.m) 



100 



Figure Al - continued 



10" 



10'- 



10" 



10" 



-■ — i i i I — 
NGC5506 



SB1. 





0.2 
0.0 
-0.2 



10 15 20 25 30 35 



10" 



10'- 



10'' 



10" 



10" 



10'- 



10'' 



10" 



10" 



SB3« 




S 


0.2 
0.0 
-0.2 



JL 



10 15 20 25 30 35 

H > > I 

SB5. 






0.2 
0.0 
-0.2 





; 1 < . J ^ 











10 15 20 25 30 35 



10" 



10'^ 



10'' 



10' 



10" 



10'' 



10'' 



10' 



SB2« 




8 

t 2 

0.2 
0.0 

-0.2 



10 15 20 25 30 35 

H , , — I 

SB4. 




10 15 20 25 30 35 



10 



Wovelength (jim) 



100 



10 



Wovelength (p.m) 



100 



Figure Al - continued 



-■ — i i i I 

ES0103-035 



SB1. 



10" 



10" 



10" 



10" 



10' 




10 15 20 25 30 35 

H , , — I 



SB3« 



10"' 



10'' 



10' 



10" 



10" 



10"^ 



10'" 



10" 



10" 



10" 



10° 




-0.2 - ■ ■ ■ : 



10 15 20 25 30 35 



H — I I I I 



H 1 1 — I I I I 



— !r- 



SB5. 





^iK^^/i^^ll i^w^^ ^ 



10 15 20 25 30 35 
_l 1 1 1 1 1 



10 



Wovelength (jim) 



100 



lO'' 



10'= 



10" 



10" 



10' 



10'- 



10'' 



10' 



10" 



10' 



SB2« 





10 15 20 25 30 35 



10 



Wovelength (jim) 



100 



Figure Al - continued 



Mrl<509 




1.4 

"> 1.2 
: o 1.0 

■ C 0.8 

: 0.6 

■ t 0.4 

■ 8:8 

0.2 t- 

r 0.0 

-0.2 



10 



15 



20 



25 



30 
j_l 



35 



10 



Wovelength (p.m) 



100 



Figure Al - continued 



10" 



10" 



10" 



10" 



10" 



"T" 



IC5063 



SB1. 




10" 



10'- 



10" 



10" 



10" 



10 15 20 25 30 35 

H , , — I 

SB3 





10" 



10'- 



10'' 



10" 



10" 



10' 



10 15 20 25 30 35 

H — I I I I 1 1 1 1 1 — I I I I 

SB5. 





10 15 20 25 30 35 



10 



Wovelength (jxm) 



100 



10" 



10'- 



10" 



10" 



10" 



10" 



10" 



10" 



10" 



10" 



SB2« 






10 15 20 25 30 35 

H ^ ^ I 

SB4. 




10 15 20 25 30 35 



10 



Wovelength (/um) 



100 



Figure Al - continued 



